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SUMMARY 


This report describes studies of CW laser -heated rocket propulsion, 
in both the flowing-core and stationary-core configurations, with most 
attention focused on the former, (It is the second such study performed 
by FSI following the work described in Ref, 1), In the present work, the 
laser radiation considered was 10, 6 [jm, and the working gas was unseeded 
hydrogen. The areas investigated included initiation of a hydrogen plasma 
capable of absorbing laser radiation, the radiation emission properties of 
hot, ionized hydrogen, the flow of hot hydrogen while absorbing and radiating, 
the heat losses from the gas and the rocket performance. The stationary- 
core configuration was investigated qualitatively and semi-quantitatively. 

It was found that the flowing -core rockets can have specific impulses 
(I S p) between 1500 and 3300 sec. They are small devices, whose heating 
zone is only a millimeters to a few centimeters long, and millimeters to 
centimeters in radius, for laser power levels varying from 10 to 5 000 kW, 
and pressure levels of 3 to 10 atm,' Meat protection of the walls is a vital 
necessity, though the fraction of laser power lost to the walls can be as low 
as 10% for larger powers, making the rockets thermally efficient, 

A number of major areas of uncertainty have been identified for 
further exploration. Chief among these is the properties of laser supported 
combustion (L.SC) waves in hydrogen. Others are an efficient method of 
heat protection, the effect of radial temperature profiles on radiation loss, 
and the effect of wider variation of operating parameters. For the stationary- 
core concept, absorption lengths of hot hydrogen radiation in cold hydrogen 
are needed, and a study of entrainment of the core gas by the flowing 
propellant gas should be made. 
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1* INTRODUCTION 


This is the final report on the second contract issued by NASA/ 

Lewis Research Center to Physical Sciences Inc. for a study of a laser- 
heated rocket. The results of the first contract were reported in Ref* 1. 

The present study dealt with the details of a number of aspects of a 
laser-powered rocket engine absorbing CO^ laser radiation (10,6 \im) by 
inverse Bremsstrahlung, and using hydrogen as the working gas. There 
were two general areas of study. The first was the initiation of a high 
temperature plasma capable of absorbing the laser radiation* The second 
was the analysis of a steady coaxial-flow, laser-powered rocket engine in 
two different configurations* One configuration was a flowing core of plasma, 
absorbing the laser energy, being heated, and expanding to generate thrust. 
This core is surrounded by a cooling buffer gas flowing the annulus between 
the core and the engine walls. The other configuration was a stationary 
stable plasma absorbing the laser energy, and transmitting it to a propel- 
lant gas flowing around it, which in turn expands tc generate thrust. 

The general geometries of these two configurations are shown in 
Fig. 1, 1, In the flowing-core concept, cold hydrogen flows into the engine, 
parallel to the laser beam, and flows into a standing laser supported 
combustion (LSC) wave. This wave heats the hydrogen up to a temperature 
at which it can absorb laser energy by inverse B r em s st r ahlung . The hydro- 
gen is further heated by this process while flowing in a variable area 
streamtnbe, until all the laser energy is absorbed* The temperature of the 
hydrogen is sufficiently high that it is a dissociated and ionized plasma behind 
the LSC wave, and radiates energy. When all the laser energy is absorbed, 
the heated plasma expands through a nozzle to generate thrust* This 
plasma core streamtnbe will transmit energy through its side walls by both 
convection and radiation, and to protect the rocket walls a eo- flowing buffer 
gas stream not heated by the laser may be necessary. The amount of buffer 
gas flow necessary is determined by the amount of energy lost radially by 
the plasma core. 

In the stationary- cor e concept, a hot core of stationary hydrogen 
continually absorbs laser radiation and re-radiates it to a surrounding 
annular streamtube of flowing hydrogen propellant gas* The propellant gas 
absorbs the core radiation and becomes heated* After passing the core, it 
is expanded to generate thrust. The propellant gas also acts as a buffer to 
protect the walls from the hot core. 
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FLOWING -CORE CONFIGURATION 




STATIONARY - CORE CONFIGURATION 


Fig. 1.1 Geometry of Laser-Heated Rockets 
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Clearly, these kinds of engines depend on a large number of physical 
processes. There is laser radiation absorption, for which the Inverse 
Bremsstrahlung mechanism has already been defined as the desirable 
mechanism in Ref* 1, There is the radiation emission by the hydrogen 
plasma, which is a loss in the flowing-core concept and the main energy 
transfer mechanism in the stationary- core concept. There is the absorption 
of plasma radiation by colder hydrogen, which heats the propellant gas in 
the stationary-core concept, and the on-coming gas in the LSC wave. This 
absorption also can protect the walls. There is convective heat transfer 
from the hydrogen plasma core to the surrounding environment. There is 
the flow of heated hydrogen in a variable area streamtube while dissociating 
and ionizing. There is the LSC wave mechanism by which a laser can 
convert room-temperature hydrogen to a plasma which can absorb the laser 
energy* There is the interfaeial mixing as the propellant gas flows past the 
core in the stationary- core concept. 

In addition to these processes important to the steady state operation 
of the rocket, other physical phenomena are pertinent to the initiation of the 
absorbing plasma. These include the possible mechanisms of gas breakdown 
by electric discharge, by a focussed laser, or by absorption of laser energy 
by molecules or particles. 


Although all these processes are involved in the operation of a laser- 
heated rocket, not all of them are studied in detail in this report. Here, the 
emphasis is on the flowing-core concept, and the processes important to 
its operation are the ones to which most attention is given. The work can 
be viewed as an extension of the studies presented in Ref. 1. 

Information is presented on the three mechanisms of plasma initia- 
tion by electric arc discharge, by laser-induced breakdown, and by molecu- 
lar and particulate abosrption, in Secs, 2, 3 and 4. A detailed study of the 
radiative emission from hydrogen is given in See. 5, leading to expressions 
useful for coupling the emission to the fluid mechanics of flowing hydrogen. 
Section 6 presents an estimate of the intensity vs wave speed for LSC waves 
in hydrogen by a rough extension of one of the theories for air. In See, 7 
is a detiled description of a hydrogen plasma in chemical equilibrium, 
absorbing laser radiation by inver se Bremsstrahlung, radiating, and flow- 
ing through a variable area channel. Section 8 presents calculations on the 
heat loss by both radiation and convection of the flowing plasma core, and 
an estimate of the radiation which can be absorbed by particles in a buffer 
gas. In. See. 9, the performance of the flowing plasma core as a rocket is 
given, based on the description of the core in See. 7, The stationary-core 
rocket is discussed in Sec. 10 in a qualitative and semi-quantitative way, 
based on estimates of the physical phenomena involved. The conclusions 
obtained so far in the study of laser-heated rockets are presented in Sec, 11. 
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The cases calculated for the flowing -core streamtube were laser 
powers of 10, 100 and 5000 kW, each at 3 and 10 atm initial pressure. 

These choices were made to provide information on sizes of engines in which 
NASA/ Lewis has expressed interest. An experimental engine designed to 
operate at 10 kW, 3 atm is being designed and built by Rocketdyne Division 
of Rockwell International, under contract to NASA/Lewis, and the same 
organization is doing a study of a 5000 kW engine* PSI has aided Rocketdyne 
with the fluid mechanics and radiation aspects of their program, under a 
subcontract. The choice of a 100 kW laser was made to provide an inter- 
mediate size between 10 and 5000 kW, which might some day be tested. The 
calculations at the higher pressure of 10 atm were made to show the effect 
of pressure level on the core flow. 


In addition to the laser power, and pressure level, there are a 
number of other parameters which could be varied in the calculations. 

One is the laser intensity which determines the mass flow rate per unit 
area into the core. In the present work this was held fixed at 3. 67E 9 W/m 
(3. 67E 5 W/cm^), which see rms a value at which the LSC wave could be 
supported. When better knowledge of JLSC waves in hydrogen is available, 
other values may be of interest* (The cross-sectional area of the rocket 
is inversely proportional to this intensity for a fixed laser power. ) A 
second parameter which could be varied is the streamtube velocity distribu- 
tion, which was also held fixed in the present calculations. Such parametric 
calculations could be performed by a slight improvement of the core flow 
program used here, and would aide in mapping out the regions of interest 
for laser -heated rockets. 


2. PLASMA INITIATION BY ELECTRIC DISCHARGE 


One possible technique for initiating laser supported combustion 
(LSC) waves in hydrogen gas is through use of an electric discharge- The 
purpose of the discharge is to supply an initial density of electrons in the 
gas so that the degree of absorption of laser energy by inverse Bremsstrah- 
lung is sufficient to drive the gas to the required temperature for formation 
of a stable LSC wave* Electric discharges have previously been success- 
fully used to initiated LSC waves in air^, No such studies are available 
for hydrogen- Thus the first question to be answered is what type of 
discharge is best suited for LSC ignition* 

Electrical discharges can be uroken into two basic groups, depending 
on whether they are or are not self-sustaining. In the latter group, desig- 
nated as Townsend discharges, electrons are formed at t ie cathode by a 
technique such as thermionic emission,; these are then accelerated by the 
electric field and, depending upon the field strength, these initial elections 
can further ionize the gas within the discharge gap. The characteristic 
electron densities in such discharges are ~ 10^/ cc and the electron temp- 
eratures can be several electron volts, well out of equilibrium with the 
gas translational temperature. 

As the voltage is increased across the electrodes a point is reached 
where a spark is created across the gap. This spark is a result of gas 
breakdown and will occur even if the cathode does not emit electrons, thus 
the term self-sustaining. The voltage at which a spark will occur, V s , is a 
function of the product of gas pressure and gap distance, pD. The functional 
dependence of V s on pD is described by a Pasehen curve which is peculiar 
to the discharged gas. The Pasehen curve for hydrogen is shown in Fig. 2* 1 
and, as an example, it can be seen at 1 atm and a gap length of 1 cm an 
applied field over 10,000 volts is required for sparking, 

There are a variety of types of self- sustained (spark) discharges^; 
however, only two characteristic forms will be considered in the present 
discussion. The first of these is a glow discharge, characterized by elec- 
tron densities of ~ 1015/ ec and electron temperatures of several eV, and 
the second is a thermalized arc discharge characterized by electron 
densities of ~ 10^/ee and an electron temperature of order one eV, in 
equilibrium with the gas translational temperature. 
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These three types of discharges, Townsend, glow and fchermalized 
arc, have been examined under generic conditions for discharge initiation 
in order to determine which type is best suited for LSC ignition. This is 
not meant to imply that all three types of discharge may be attainable in a 
specific environment* Indeed Townsend and glow discharges are more 
favorably produced at sub" atmospheric pressures; however, the properties 
of a discharge are determined not only by pressure but also by numerous 
other factors such as voltage, gap size, electrode configuration., and the 
energy content of the relevant electronics. 

The scenario considered is that of a discharge operating perpendic- 
ular to the gas flow, the latter being eo- axial with the laser beam. The 
distance between anode and cathode is not specified but must be larger than 
the laser spot size since laser irradiation of the electrode elements is not 
desirable. The discharge length and gas velocity have arbitrarily been 
chosen to be 1 cm and 10^ cm/ sec respectively, so the gas flow time in 
the discharge region is 1 msec. Note this assumption does not imply that 
the discharge occurs simultaneously along the 1 cm length. In the ease of 
a spark discharge the streamers would be considerably smaller than this. 
However, once a spark is initiated, it will attempt to travel along the path 
of least resistance and thus follow the flowing gas for some distance. 

Once the discharge is initiated the gas in the discharge region will 
begin to absorb the laser energy via neutral- electron inverse Iremsstrah- 
lung (or if the electron density is high enough by ion- electron inverse 
B r em s s t ra hlung ) . The total energy absorbed per unit volume of ionized 
gas is given by 


A E 


I k L I * 


( 2 . 1 ) 


where I is laser intensity and k ^ is the inverse neutral Bremsstrahlung 
absorption coefficient at the laser frequency, defined by 

k L> = a L nn E cm ~ 1 (2. 2) 

where n is total gas number density, n E is the electron number density and 
a E (T) is a function of temperature (or, more appropriately, electron 
temperature in non- equilibrium situations) and frequency which is specific 
to the absorbing gas. 
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As the gas absorbs energy it will neat, causing an initial pressure 
rise and subsequent expansion which complicates the evaluation of relation- 
ship (2. 1). In the present analysis a firm upper bound on the energy 
absorbed by the gas has been determined by assuming that the number 
density of the gas remains constant as the temperature increases. In this 
ease, since the electron, density of the discharge is specified, k will be 
determined by the form of (T). In the cold gas the dominant neutral 
species will be however, as the gas is heated, dissociation will occur 
and hydrogen atoms will become the dominant neutral partner for the 
absorption process. The reduced absorption coefficients kjJ for H and M 2 
are shown in Fig. 2.2 as a function of temperature, for a wavelength of 
10,6 pm. These predictions have been developed from sources described 
in Reference 1. It is to be noted that the reduced absorption coefficients for 
the two gases are quite similar and that they scale approximately inversely 
with the temperature. The reduced absorption coefficient may be related to 
ajj (T) of Eq. (2, 2) by the relationship 


k L ' = a L (T) kT , 


(2.3) 


where k is the Boltzmann constant, and thus for these gases a^(T) is 
approximately a constant for X = 10.6 pm, Therefore, from Fig. 2.2, 
may be approximated by 


-37 -1 

k T 3. 5 x 10 tin_, Gm 

Mr 


(2.4) 


Taking n = 2. 5 x 10 part/cc, corresponding to p = 1 atm and T = 30O°K, 
Eq. (2. 1) then reduces to 


AE = 10 n I, Joules /cc 

IL 


(2.5) 


Equation (2. 5) provides an approximate estimate of the amount of 
energy absorbed in one msec by a constant density, constant electron density 
plasma. The increase in temperature of the gas which is a result of this 
absorption maybe deduced by consideration of the thermodynamic properties 
of heated hydrogen. The variation with temperature of the enthalpy /unit 
mass of constant pressure hydrogen plasmas is shown in Fig. 2. 3 for 
pressures between 1 - 1000 atm. These predictions were taken from the 
analysis of Pateh^. The enthalpy variation of a constant density hydrogen 
plasma may be deduced from a cross interpolation of the constant pressure 
data. Thus for a given amount of absorbed energy, as predicted by Eq, (2. 5), 
a unique plasma temperature may be determined. 



Fig a 2.2 Central Inverse Br ems straUong Reduced 
Absorption Coefficient for 10,6 pim 





Equation (2,5) has been evaluated as a function of laser intensity for 
the three generic discharges discussed earlier. The results are shown in 
Fig* 2.4* The final hydrogen plasma temperature which would be attained 
for a given amount of absorbed energy is shown on the right hand scale* 

Also shown as tick marks on the curves is the temperature at which the level 
of thermal ionization exceeds the ionization level created in the discharge. 
This latter temperature is a useful criterion for ignition inasmuch as a stable 
LSC wave would exhibit ionization levels higher than that created in the 
discharge (at the specified total densities). It can be seen from Fig. 2.4 
that a Townsend discharge will not produce ignition and that a glow discharge 
would be effective at best only for intensities > 10^ W/ern^, Recalling that 
the results of Fig, 2*4 represent a strong upper bound on the absorbed 
energy, it would appear that the thermalized arc discharge is best suited 
for the successful electric discharge initiation of ESC waves in hydrogen. 

From a systems viewpoint it is desirous to define the minimum 
discharge requirements for successful ignition, in terms of electrical 
energy content, device size and sophistication. Unfortunately, it is difficult 
to provide such a definition by theoretical means. The production of stable 
thermalized arcs is still more of an art than a science. One standard tech- 
nique for their generation is through use of an over- voltaged capacitor , The 
evidence 5 suggests that a plasma with electron densities of ~ 10^ 5 /ee is 
forced from the cathode end by high velocity ionizing waves into a spark 
channel* This streamer provides a finite resistance path for discharge of 
the capacitor. As the capacitor discharges, the gas in the streamer region 
will be heated by ohmic power dissipation. Presently, this behavior cannot 
be predicted a priori inasmuch as the plasma resistivity is a function of 
plasma electron density and streamer size, and the electron density in turn 
is connected to the rate of ohmic heating. This analysis is further compli- 
cated by the hydrodynamic behavior of the gas**. The initial spark can lead 
to formation of a cylindrical shock wave which can produce further ionization. 

In any event, such discharges are readily produced in hydrogen and 
thermal equilibration times have been observed to be tens of nanoseconds 5 * 
The discharged gas is low density (relative to p*) but highly ionized. Indeed 
measurements^ performed in hydrogen in near atmospheric discharges 
exhibit electron concentrations near the anode as high as 1 0 ^ ® /e e and final 
temperatures as high as 50 - 60, 000°K. It is concluded that LSC wave 
ignition may be achieved through use of a thermalized arc discharge; how- 
ever, it is felt that optimum discharge conditions should be determined by 
laboratory experiments rather than through theoretical analysis. 
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Fig. 2.4 Laser Energy Deposition for Various Discharge 3 

Electron Densities (Gas Density 4.9 x lO 1 ^ parfeicles/cm ) 
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3. PLASMA INITIATION BY LASER-INDUCED BREAKDOWN 


3. 1 INTRODUCTION 

Free electrons can absorb laser radiation during collisions with 
heavy particles such as ions, atoms or molecules* In this process, known 
as inverse Br emsstrahlung, the role of the heavy particle is to conserve 
momentum. Little of the energy of the absorbed photon is transferred to 
the heavy particle. Free electrons transfer energy among themselves 
easily during collisions; therefore the laser energy absorbed by the elec- 
trons is rapidly distributed among all the free electrons. The electron 
distribution can then be characterized by an electron energy 0 roughly equal 
to two-thirds the average electron energy just as kT is related to two-thirds 
the average energy in a Maxwellian distribution, (The factor of two-thirds 
relating e and the average energy is exact only for Maxwellian distribution; 
nevertheless the characteristic energy e is a useful concept* ) 

Energy transfer from the electrons to the heavy particles is less 
efficient than electron-electron transfer; hence it is possible for the gas 
temperature to lag behind the electron energy. For strong fields, electron 
heating occurs so rapidly that the few free electrons reach high characteris- 
tic energies, that is, a sizable fraction of the heavy particle ionization 
potential, without heating the gas noticably. At such energies a significant 
number of the electron - heavy particle collisions result in ionization. A 
chain reaction then occurs with both the original electron and the newly 
created electron absorbing laser energy and participating in further ionizing 
interactions. After a short time the avalanche created by electron - impact 
ionization leads to a sizable concentration of high energy electrons in a 
otherwise cold gas, i. e. , to gas breakdown. 

Experimentally, as laser intensity is increased from an initial low 
level an abrupt threshold is passed above which this breakdown occurs. ® 
This threshold intensity can be estimated by balancing two important physi- 
cal mechanisms. The first balance determines the number of electrons 
which remain in the laser - heated region. The rate of production of new 
electrons must equal the rate at which they are lost* The second balance 
determines the characteristic energy the electrons reach. They are heated 
by the laser, but lose heat by collisions with the gas molecules. From an 
analysis of there two balances, an estimate of the intensity for laser break- 
down in hydrogen will be obtained. 
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3. 2 ELECTRON PRODUCTION AND LOSS 

In addition to production of electrons through electron - impact 
ionization, electrons may be freed by multiquanta absorption, the process 
in which an atom or molecule absorbs many low energy photons rather than 
one high energy photon. Whenever an appreciable number of electrons is 
present electron - impact ionization dominates. Multiquanta absorption is 
important only as a potential source of the initial electron population 
required for electron - impact ionization to dominate. ^ 

In molecular gases the loss mechanisms generally are attachment of 
an electron to a neutral molecule, recombination of an electron with an 
ionized species (with or without the presence of a third particle) and electron 
diffusion out of the volume influenced by the laser. Attachment is unimpor- 
tant in hydrogen. The recombination rate is also small, particularly during 
the initial stages of breakdown. Diffusion is the predominant loss process. 

The net rate of electron production in hydrogen can thus be expressed 
as 


s v.n E ♦ r • dP» e . (3.1) 

where m-jj* is the number density of free electrons, is the ionization rate 
per electron (the concentration of neutral molecules is absorbed into the 
definition of v.) and D is the diffusion coefficient. 

For rough estimates the diffusion loss mechanism can be approxima- 
ted by 


Dn. 


V D n E = V ' (S ^ V 


E 


R 


b 


(3.2) 


where R^ is the radius of the laser beam. As long as the electron density 
is low the electron diffuses freely. At high concentrations diffusion proceeds 
at the slower ambipolar rate. For our estimates of threshold intensity we 
will use the free diffusion eoeffient. 

Breakdown requires electron production to exceed diffusion losses. 
The threshold electron energy for breakdown can be estimated by balancing 
losses and gains. 
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(3.3) 


V i * 6) = V D (6) 


The quantities Dn^ and e, where nj^ is the number density of free 

molecules, are known aB functions of the electric field per particle, E/nj^j, 
from experiments with static electric fields. U, ^ Despite the differ- 
enee in the electron distribution function for DC fields as compared to AC 
fields calculations have indicated that the ionization rate for the same € is 
very similar* We have thus used the DC experimental radius to determine 
v i^ n M an< ^ ^ n M as functions of electron energy e. ^ In Fig. 3. 1 we have 
plotted versus e the diffusion frequency per molecule for an electron, 
v D/ n M f ^ or var i° us values of pR^ ( (where p = (3. 034E-17) nj^ at 293K), 
and the ionization frequency per molecule for an electron, The 

extremely steep increase of v^/nj^ with e indicates why the threshold is 
abrupt. A small increase in e above the temperature needed to balance 
diffusion losses results in ionization completely dominating diffusion* For 
example ionization and diffusion balance at e ~ 2* 7 eV when pR^ = 1 torr- 
cm. If the electron temperature reaches 2, 8 eV for the same pR^, we find 
that 


v. 
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v 
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cm 


/sec 


(3.4) 
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For 1 atmosphere hydrogen at 293K, ** 2. 5 x 10 cm so that solving 

the electron production Eq. (3* 1) leads to 


n 


n 

o 


exp (5.0 x 10 



(3b 5) 


Thus if the breakdown requirements are exceeded by even a small amount 
the gas breaks down rapidly. [Of course other time scales are also 
involved - the time to heat the electrons and the time required to produce a 
reasonable initial electron population if one is not present.] 

3. 3 ELECTRON ENERGY BALANCE 

The minimum laser intensity necessary to reach the breakeven 
election energy is determined by equating the electron heating by the laser 
and the power lost through collisions with the molecules. The electron 
heating is described by 
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Fig. 3,1 Diffusion Frequency and Ionization Frequency 
in Hydrogen at Various Values of 
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where e and are the charge and mass of the electron, respectively, (u 
is the angular frequency of the laser radiation, E is the RMS electric field, 
v m momen i :urn transfer collision frequency per electron and Vi n the 
inelastic collision frequency per electron. The first expression on the 
right hand side represents inverse Bremsstrahlung which can be calculated 
from collision theory. For our purposes, however, it is simpler to use 
the expression given above and the experimentally determined values of 
v m and Vi n versus e as determined from DC experiments. 

The minimum laser intensity 1 necessary to reach e is found by 
setting the right hand side of Eq. (3. 6) equal to zero. 


We have 


I = cE /4ft 


\ n M / V m ' 
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(3.7) 


For I in watts /cm and e in eV, Eq. (3.7) becomes 


, . . . 18 ^ in^M* , 2 . 2 V „ 

I " 1,5xl ° (v /n..) (U) +V m )S * 

m M 


(3.8) 


It is now possible to investigate qualitatively how close the breakeven inten- 
sity is to the true threshold intensity. In the eV range v^ n increases 
extremely rapidly with small increases in e where as v m is almost constant. 
Hence the second term of Eq# (3*6) is negligible throughout most of the 
heating. It serves only to stop the heating abruptly at the final electron 
energy. Because of the sensitivity of both the production rate and the heat- 
ing rate to small changes in e, the two balances, Eqs* (3. 3) and (3*6) 
should lead to accurate estimates of the threshold intensity. There is 
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greater uncertainty in the experimental values of v , Vj, Vi n , Dnj^ and 
e then in the identification of the threshold intensity with the breakeven 
intensity. 


3.4 BREAKDOWN ESTIMATE 

For a given beam radius R^ and pressure p, the laser intensity to 
produce breakdown can now be estimated. Using Fig. 3. 1, the energy e 
at the crossing point of the Vj7n^ and Vp/nj^ curves is found. Using this 
energy, the quantities v^ n and V m are obtained from Refs. 10 or 11* Then 
the laser intensity I is calculated for the desired circular frequency m 
from (3. 8). Using this procedure the threshold intensity necessary to 
achieve breakdown in pure hydrogen initially at 293° K and 3 atmospheres 
pressure is plotted versus laser beam radius for a 10.6 laser in Fig* 3.2. 

The threshold intensity is seen to decrease with increasing spot 
size until a plateau near I - 4,5 x 10^ W/em^ is approached. The decrease 
is caused by the lowering of the diffusion loss rate with increasing radius. 
The plateau represents the minimum intensity needed to raise the electron 
energy high enough for appreciable ionization to occur. The intensity of a 
10 kW laser is also plotted as a function of beam radius. Over a limited 
range of pressure and wavelength the breakdown curve scales as follows: 
for a given laser intensity and wavelength, the beam radius is inversely 
proportional to pressure; for a given radius and pressure the threshold 
laser intensity is proportional to the wavelength squared. 

Focussing a powerful laser unto a small volume will produce a high 
concentration of free electrons within a small volume. The effectiveness of 
laser induced breakdown in initiating an ESC wave is probably not limited by 
the electron concentration achievable but rather by the volume over which 
such a concentration can be sustained. 

3. 5 LASER-INDUCED BREAKDOWN WITH PARTICULATES 

The effect of particle contaminants on the breakdown threshold for 
air is a subject which has recently received considerable attention. It has 
been found by numerous investigator s 14 - 16 that the intensity threshold in 
the presence of particulates may be lowered as much as three orders of 
magnitude from the clean air threshold. Although there is no data for the 
effects of particulates upon the breakdown threshold in hydrogen, the mech- 
anisms by which partieules induce breakdown at lower intensities in air will 
also lead to lowering of the threshold in hydrogen* 
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Fig, 3.2 Laser Intensity Required for Breakdown in Hydrogen 
and Laser .Intensity Available from 10 kW Laser 
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A summary of the observed breakdown data in air is presented in 
Fig, 3,3. The air breakdown threshold is a function of pulse time, beam 
diameter and particle size. The lowest observed threshold was of the 
order of a megawatt per em^; however, 30 \lxn particles were required to 
initiate the breakdown. In addition, it was observed^® that no breakdown 
occurred with the larger particles at the lowest intensities until lO^isee 
into the pulse. Shorter pulses do not yield breakdown with contaminants 
unless the flux is considerably higher. 

The physics of particle induced breakdown is not well understood. 

For short pulses the particles evaporate rapidly, and cascade breakdown 
occurs in either the particle vapor or the air surrounding the particle. 
Triplett and Boni*® have performed a theoretical study of the interaction 
of laser radiation with a single suspended particle. Their analysis applies 
only when the laser intensity is sufficiently large that the rapidly evapora- 
ting or "exploding" particle drives a strong shock wave into the surrounding 
air. Breakdown then occurs via a non equilibrium cascade process (des- 
cribed above) in either the vapor or shock heated air. Calculations indicate 
that reductions in the breakdown threshold to slightly greater than one order 
of magnitude below that of clean air are feasible by this mechanism. It is 
not possible to predict with this mechanism the reduction in the threshold by 
the observed several orders of magnitude at longer laser pulses. In addi- 
tion, an "incubation" time or delay time between laser turn-on and break- 
down was observed with the long laser pulses^ ( 10 ^isee). Since the delay 

times are 10-30 ^lsee, the breakdown mechanism is nonequilibrium cascade 
ionization of shock heated air or vapor. It appears that the vapor formed 
by the evaporating particle heats in local thermodynamic equilibrium and 
transfers its energy to the surrounding air by conduction and radiative 
transport. The air ultimately heat via inverse Brems strahlung absorption 
until the formed plasma becomes opaque to laser radiation. However, at 
present there are not theoretical models capable of predicting the break- 
down threshold under these conditions. 

Therefore, although, it is not possible to predict the breakdown 
threshold in particle seeded hydrogen, we would expect this to be a viable 
technique for igniting the plasma at lower intensities than predicted by 
Fig, 3,3. 
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4. 0 PLASMA INITIATION BY MOLECULAR AND 
PARTICULATE ABSORPTION 


The two methods of initiation discussed in the previous two sections, 
electric discharge and laser-induced breakdown, provide the electrons to 
initiate inverse B r em s s t r ah lung absorption by a nonequilibrium process. 

A third alternative is to provide the electrons in an equilibrium process by 
heating hydrogen into its ionization region. For pure hydrogen, that means 
heating to 10,000 K. If the hydrogen is seeded with an easily ionizable 
material, a much lower temperature would suffice. Cesium, for example, 
ionizes around 3000 K. 

Such heating could be accomplished if one could find molecules or 
particles which would absorb 10. 6 pm laser radiation, become hot and 
transfer their heat to the gas. In Ref. 1, Sec, 3, a study was made of 
possible molecular absorbers for 10. 6pm radiation. It was concluded that 
high temperatures could not be reached because the molecular absorbers 
dissociate (and so stop absorbing) at a few thousand degrees. No mole- 
cules which circumvent this limitation were identified in Ref. 1, nor have 
any come to light during performance of the work reported here. However, 
if seeded hydrogen is considered, then the gas need only be heated to a few 
thousand degrees for ionization to occur. Molecular absorption might be 
an effective means of accomplishing this heating. 

Particulate absorption was also studied in Ref. 1, Sec. 4, af well 
as in Ref. 19. Here the particles absorb laser photons, get hot, and con- 
duct heat to the hydrogen. Here again a temperature limitation was found. 
An upper bound for particle vaporization is 5000 K, limiting the temperature 
which can be reached to below the pure hydrogen ionization region, but 
within the ionization region of seed materials. 

There was noted in Refs, 1 and 19 some possible benefits from hot 
particles, other than heating the gas. The vapor produced from hot parti- 
cles may itself absorb the laser radiation, and may also have a lower 
ionization potential than the hydrogen. Furthermore, small particles can 
emit electrons thermionieally as they are heated, though this emission is 
inhibited by electrostatic effects as the particle charge builds up. However, 
the electrons produced may be sufficient to support laser-induced break- 
down. 


Nevertheless, it seems at the present time that molecular or parti- 
culate absorption is not a favorable scheme for initiating an un seeded 
hydrogen plasma, compared to electric arc or laser-induced breakdown, 
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5. RADIATIVE EMISSION FROM HYDROGEN 


5. I INTRODUCTION 

As the laser heats the hydrogen to high temperatures the radiation 
from the gas increases dramatically* In an ionized plasma, radiation con- 
tributes importantly as an energy transfer mechanism within the plasma 
and as a loss mechanism from the plasma. The radiative transfer within 
the gas affects both the axial and radial temperature gradients. Thus a 
precise description of the physical state of the gas in the rocket core 
requires that the radiation and fluid mechanics be coupled. It is beyond the 
scope of this report to deal with all aspects of the coupled problem. In 
keeping with approximations to be made in the fluid mechanics, such as 
the omission of the radial temperature gradients, we neglect radial radia- 
tive transfer within the plasma. Further-more, the axial radiative trans- 
fer within the hot gas is expected to be small compared to the laser heating; 
it will not be considered here. Axial radiative transport from the hot 
plasma to the cold incoming gas may be important* Its effect is included in 
the LSC wave properties discussed in the next section* 

The only plasma-produced radiation remaining to be considered is 
the radial flux leaving the gas. This radiative transport enters into 
rocket engine analysis in two ways. First, the power lost radially must 
be absorbed somewhere else in the rocket, either in the walls or in a buffer 
gas. Second, the radial radiation which escapes the plasma must be 
accounted for in the fluid mechanical equations as a loss. At sufficiently 
high temperatures, in fact, it is the major loss mechanism and effectively 
limits the maximum temperature attainable. 

The approach we follow here is first to study the local radiation 
properties of the plasma, and then to develop, using the local properties 
as guides, a reasonable approximation to the coupled radiating gas dynamic 
system. 


5, 2 LOCAL SPECTRAL RADIATION PROPERTIES 

The spectral radiation properties of a hydrogen plasma are deter- 
mined by the local state of the system. We assume that the plasma is in 
local thermodynamic equilibrium, so that temperature and pressure define 
the state of the gas. The concentrations of the various species is determined 
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through the equilibrium conditions* The temperatures and pressures of 
interest are those expected in a rocket heated by laser absorption; that is, 
temperatures between 10, 000 K and 20, 000 K and pressures between one 
and ten atmospheres. In the interest of simplicity we need only include 
effect' which are expected to contribute 10% or more of the total radiation. 

The basic ingredients for determining local radiative properties are 
the cross sections for the possible interactions and the number density of 
the various species. The use of local thermodynamic equilibrium assumes 
that collisions among the various species of the gas are frequent enough 
to ensure an equilibrium distribution of velocity and of species, corres- 
ponding to the local temperature and pressure. With this assumption and 
the principle of detailed balance the radiative transfer equation can be 
written in a simplified form involving only one function depending upon the 
concentration of species and the interaction cross-section. The radiative 
transfer equation has the following form:^ 


^ • V \ (T, p) [B^ (T) - I (x, n) ] (5. 1) 

where ui - hp v is the photon energy; 1^ is the spectral radiation intensity 
which depends on u>, x, and the direction specified by the unit vector n; 
and (T) is the equilibrium spectral radiation intensity given by 


B (T) = 
(!) 


U>... 


2 ~ 3 , UJ 
e n (e 



(5.2) 


where c is the velocity of light, h is the Planck constant divided by 2TT and 
k is the Boltzmann constant. The effective absorption coefficient is related 
to k (T,p), the absoption coefficient uncorrected for reabsorption, by 


» -w/kT 

k^ (T, p) = (1 - e ) k y (T, p) . {5. 3) 

The absorption coefficient k (JU (T, p) is, in turn, easily expressed as a sum 
of products of absorption cross sections and species concentrations. The 
temperature and pressure dependence of k^ enter mainly through the con- 
centration of species. However , secondary temperature variations occur 
whenever the cross section must be averaged over the relative flux of some 
of the involved species, as occurs in inverse B rems strahlung . 
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The first step in determining k fjj (T,p) is to detail which interactions 
are expected to be important for the range of temperatures and pressures 
expected. The interactions to be considered, in turn, depend upon the 
species which are abundantly present. When radiation is an important loss 
mechanism the concentration of molecules is several orders of magnitude 
lower than the atomic concentrations. We therefore consider a plasma in 
which only atoms, excited atoms, ions (protons) and electrons are present. 
Their relative abundance is determined through the Saha equation 


n E n I 


n, 


Q el I Q el E 


el A 



(5. 4 ) 


where n-g, nj and are respectively the concent ration of electrons, ions 
and hydrogen atoms, QelE* Q e li Q e lA aire partition 

functions of the electron, ion and the neutral atom and 0j is the ionization 
portntial expressed in temperature units, (In more accurate calculations 
than are necessary here, ©j would be replaced by - A ©j where A ©j is 
the effective lowering of the ionization potential corresponding to termina- 
ting the sum in the electronic partition function of the neutral atom after a 
finite number of excited states). 

The transitions which lead to emission can now be listed# They are: 

1) capture of a free electron by an ion 

2) B r em s st r ahlung involving electron* ion collisions 

3) IB r em s s t r ah lung involving electron-atom collisions 

4) radiative decay of an excited atom# 

5. 3 CONTINUUM TRANSITIONS 

The spectral properties of the first three types of transitions, the 
transitions involving the electron continuum, can be estimated by semi- 
classical formulas. The quantum corrections to these formulas need only 
be included when the quantum correction is large and the radiation emitted 
is a sizeable portion of the total radiation produced. 

From Kramer f s semi- classical analysis of electron-proton Brems- 
st rah lung the absorption coefficient is given by^ 
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where e is the charge of the electron and the pressure dependence enters 
only through the number densities. 

Similar semi-classical formulas can be derived for bound-free 
absorption wherein an electron of energy E n in the nth quantum state of 
an atom absorbs a photon of energy u> resulting in a free electron kinetic 
energy uj + g and a free ion. The semi-classical cross section for absorp- 
tion is^2 

a 10 
8tt e m E 

V n = TlrTr? ^ s'" (5-6 

3 ,/ 3 ( n ) c uj n 

The total bound free absorption coefficient for radiation of given uj is then 


n <j 
n co > u 


n & n 


where n is the number density of atoms in the nth quantum state, and where 
n^"” is the lowest quantum state which can be ionized by a photon of energy a); 
that is 


❖ 2 
{n ) 
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The Bremsstrahlung due to ion- electron collisions Eq. (5. ) can be 
re-expressed in a form similar to Eq. (5. 7) by using the Saha Eq. (5.4) 
with Q e u = 1 to convert the product ng Uj into the concentration of ground 
state atoms, n^Q and by similarly expressing n n in terms of ii^q: 
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(5.9) 
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(5.10) 
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The sum of the absorption coefficients for Bremsstrahlung and free- 
bound transitions is 
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where we define x n = — 


W 


, and we have used in (5. 5) the expression for the 


n T 


ionization potential of hydrogen in terns of fundamental constants, 


9j = e 4 nig /2 k F 2 


(5. 12) 


to get the first term in brackets in (5, 11). 
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If U) is given in eV and n^Q in em , then k m in cm is 
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(if the lowering of the ionization potential were important, Eq. (5. 11) would 

A6 t /T , 

have 1/2 replaced by e A / 2x^). 


The emission coefficient (T, p) is defined as 


e (T, p) = k (T, p) B 
0) 0) 0) 
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The semi- classical approximation to (T,p) found from Eqs. (5. 11) 
and (5. 13) is shown in Fig. 5, 1 for p = 3 atm and T = 14, 000 K. The total 
continuum power emitted per unit volume, neglecting reabsorption effects, 
is found by integrating e^, over solid angle and over the photon spectrum. 
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P cc (T ’ P) =/*> J\ (T, P) % = 4tt J (T,p) d 
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(5. 14) 
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In Fig. 5*2, IP (T, p) is plotted versus temperature for a hydrogen 
plasma at 3 atm. 

If, on the other hand, the vacuum ultraviolet part of the spectrum 
(u) £ 13.6 eV) is effectively reabsorbed, the continuum power lost in the 
semi- classical approximation is 


13*6 

F cc (u) 5 13o6) = 4 *r f e (T, P) dtu 
ec J uu 

o 



which is also shown in Fig. 5. 2. 

The relative magnitude of bound-free emission and Br em s strahlung 
can now be determined. The contribution of ion- electron B r ems strahlung to 
the emitted power, as compared to the bound-free emission, is estimated by 
comparing 1/xj to 2.404 when no reabsorption occurs and to 0,404 when the 
vacuum ultraviolet is strongly reabsorbed. For the temperatures of interest 
we have 0. 07 < < 0. 14* Thus Brems strahlung is not expected to be an 

important factor unless the bound free emission is inhibited by reabsorption. 

The other continuum contribution, Brems strahlung associated with 
electron-neutral collisions, is much smaller than electron-ion Brems strah- 
lung except when atoms are at least one hundred times more abundant than 
ions. ^ Only at lew temperatures, such as 10, 000 K, and high pressures 
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Fig. 5.1 Semi- classical Approximation for the Continuum 

Spectral Emissivity of Hydrogen and the Black Body 
Limits for Two Choices of Radius 






POWER LOST / VOLUME ( W/ctn 3 



- 32 - 



i 


t 


does the neutral concentration overwhelm the ion concentration sufficiently 
to make neutral Bremsst rah lung as important as electron-ion Bremsstrah- 
lung, which itself is only a small fraction of the total radiation. Moreover, 
as can be seen from Fig, 5* 2, plasma radiation is very small under these 
conditions. For our purposes, then, the contribution of neutral Bremsstrah- 
lung can be safely omitted as an important radiation loss mechanism. 

The expressions quoted above are semi-classical in origin. To 
include the quantum effects the cross- sections must be multiplied by a f 

correction factor known as the Gaunt factor. For free-free transitions at 
the temperatures and pressures under consideration the Gaunt factor 
averaged over the range of uu is approximately 1*4. (The precise value 
varies with temperature and pressure, ) Since this is a correction to an 
already small contribution it is not included in further computations. 

The bound-free transitions also have Gaunt factors associated with 
them. The most important quantum correction, about 20%, concerns the 
transition to the ground state. Rather than using a Gaunt factor, however, 
it is more convenient in this case to begin with the correct quantum cross 
section for the region of U) where the bulk of the emission occurs. From 
Bethe and Salpeter^ we have to replace the cross- section of (5, 6) when 
n = 1 by 
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We then separate off the continuum n - 1 term of the sum in (5, 11) 
and find the absorption coefficient to be 
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(5. 17) 


The last term is, of course, only to be used when < 1. This quantum 
correction is the only one we include in the calculations in Subsection 5, 6. 
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5.4 ATOMIC TRANSITIONS 

The final radiation process to be considered i6 the transition between 
two atomic bound states. The total power emitted in a spectral line corres- 
ponding to a transition from an upper state of quantum number n' to a lower 
state designated by quantum number n is easily computed from the well 
known f nn ' number for atomic transitions. 

The f numbers of the most important lines 6 are given in Table 5, 1, 
The f number of the higher lines (n 1 > 8) in the first three series are given 
by 

3 

(5.18) 

where Aj = 1.6 for the Lyman series, (n = 1), A^ = 3.7 for the Balmer 
series (n = 2), and Aj = 5.3 for the Pasehen series (n = 3). The f nn ' 
numbers for transitions involving n > 3 can be approximated as 

I = 1.96 m (~T > 

*■ \n -n' 2 

From the i * number the cross section integrated over photon energy is 
found to be 
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and the power emitted per unit solid angle is given by 
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where UJji - k ©j 
be rewritten as; 
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TABLE 5. 1 







f numbers for hydrogen for transitions 

FROM LOWER STATE n TO UPPER STATE n' 
FROM REFERENCE 25 


n 


n' 


•mm' 


I 


3 


2 

3 

4 

5 

6 
< 

8 

3 

4 

5 

6 

7 

8 

4 

5 

6 

7 

8 


.4162 
.0791 
. 0290 
.0139 
.0078 
. 0048 
.0032 
.641 
. 1 195 
.0444 
.02235 
.0124 
.0082 
.8413 
.150 
.0554 
.0271 
.0162 



The power emitted (4 TT times Eq. (5. 21)) in several of the low-lying 
lines dwarfs the total continuum contribution. In Table 5.2 we illustrate 
this feature for hydrogen at 14, 000 K and 3 atm* However reabsorption, 
which will be discussed in the next subsection, reduces the magnitude of 
the losses due to line radiation, In order to adequately analyze line trans- 
port it is essential to know not only the total energy radiated but also the 
detailed profiles of the strong lines. The profile is modified by the 
influence of other particles. In the presence of atoms, ions and electrons 
the shape of the line is broadened from the narrow Lorentz characteristic 
of the natural line shape; however the change in total power emitted is 
negligible. In an ionized plasma such as in the hot core of the rocket the 
broadening is caused mainly by the charged particles and is commonly 
termed Stark broadening* Hydrogen lines are affected more strongly by 
Stark broadening than lines of other atoms, so that in ionized hydrogen 
other broadening mechanisms such as Doppler broadening are unimportant. 

The profiles of hydrogen lines subject to Stark broadening serve as 
important diagnostic tools; as a result they have been thoroughly studied. 
Elaborate calculations of the profiles have been performed and the profiles 
carefully tabulated. ^ For our purposes, however, only the far wings of 
the Hyman lines are important. In the asymptotic far wings the profile S 
is adequately described by 2 times the Holts mark profile: 


s (a) = 2 c r /icti 
nn r 


5/2 


( 5 * 22 ) 


where C nn r is a constant and a is the reduced wavelength defined in terms 
of the line center energy a) Q , the line center wavelength X Q and the normal 
field strength 
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TABLE 5.2 

POWER EMITTED BY LINES IN TRANSPARENT APPROXIMATION. 
HYDROGEN AT 3 ATMOSPHERES AND T = 14000 K 


Power Emitted 


Line 

n> 

n 

f watts / cm.3 ”| 

La 

2 

i 

7.0 x 10 5 

14 

3 

i 

4.62 x 10 4 

L y 

4 

i 

2.4 x 10 4 

L S 

5 

i 

4.59 x 10 3 


6 

i 

2337 


7 

i 

1350 


8 

i 

866 


9 ^ « 

i 

2463 

H a 

3 

' 2 

5705' 

H|3 

4 

2 

1507 

Hy 

5 

2 

6l6 


6 

2 

314 


4 

3 

418 


5 

3 

176 


6 

3 

94.3 


5 o& 

_ _ 

80.4 

Continuum 



2.6 x 10 4 
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The profile is normalized by requiring 


fs ( a ) da = 1 


(5.23) 


Thus in our spectral plots of the emission coefficient vs u) we have 



(5.24) 


(5.25) 


The factor of two in Eq. (5,22) is required in order to include the electron 
broadening effects as well as the ion broadening. The values of C^i for 
the first three Lyman lines ^ 7 are given in Table 5, 3. Typical line 
broadened profiles are illustrated in Fig. 5. 3, Only the first four lines 
are shown isolated. The higher-lying lineB are indistinguishable. They 
merge to form an extension of the bound -free continuum. When necessary 
more detailed line profiles can be found in Griem. 26 

5.5 SIMPLE MODEL OF RADIATIVE LOSSES 

Using the spectral effective absorption coefficients found in the last 
subsection we can now proceed to model the radiation from the hot plasma 
core. 

In principle we wish to find the spectral intensity 1^ (x, n) at each 
point x and direction 'n' in the core by solving the equation, 

'n • V I (x, n v ) = k (T, p) [B (T) - I (x, e) ] (5.26) 

O 01 o o 

subject to the boundary condition that at the outer surface of the plasma the 
radiant intensity 1^ directed into the plasma be zero. Rather than attempt 
such an ambitious task here, we are satisfied to develop a model in which 
the plasma radiation can be treated as a local loss term in the gas dynamic 
equations and in which the total flux of radiation at the edges can be estimated. 



TABLE 5.3 


COEFFICIENTS FOR ASYMPTOTIC 



1 2 3. 37 x 10 

1 3 1.80 x 10 

1 4 3. 99 x 10 



m 


* 


ft 



U) (eV) 


Fig. 5,3 Stark Broadened Profiles of 1 
and the Black Body Limits 
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r 


To this end we consider the radiation lost per unit volume at temperature 
T and pressure p from the plasma to he equivalent to the average radia- 
tion loss per unit volume from an infinitely long cylinder of uniform 
temperature T, pressure p and radius r. In this way the major qualita- 
tive effects of reabsorption are estimated without including the effects of 
temperature variations. 

The power emitted per unit length at a given ts is 


P = 2if B (T) R F O' R] 
Uj (U 

where the function F is defined as 


(5.27) 


fr/2 

r( k w ‘* )=?/ d 0 cos 2 e 



(5.28) 


The function Ij and Lj are respectively the modified Bessel function of 
order one and the modified Struve function of order one. The function 
F (k^ R) gives the ratio of actual power emitted to the power emitted by 
a blaekbody of the same radius. It is plotted in Fig. 5.4. 

Two extreme limits are important. When is sufficiently small 
so that k R « 1 is satisfied, reabsorption effects are unimportant. 

The total emission per unit length is then equivalent to the emission per 
unit volume multiplied by the area; 


? £ (volume) = r2 B „, (5.29) 

I t 

The other extreme is large k . When k R » 1 reabsorption dominates 
and the cylinder looks like a blaekbody surface radiator: 


P 

UD 


(surface) 


= (ttB ) (2 ft R) = 

UJ 


2tt 2 B R 
OJ 


(5. 30) 


For comparison these two extreme limits are also plotted in Fig. 5.4 with 
the factor of Zff^RB divided out. Calculations using Eq. (5.25) to 
weight the various parts of the spectrum can be performed numerically; 
however, first estimates of the effect of radiation can be made by resorting 
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Zf- 


Black Approximation 



Fig. 5*4 Ratio of Actual Power Lost to Black Body Limit on Power 
Lost for an Infinite Cylinder* The Transparent and Black 
Body Approximations Are Also Shown. 


8 . 

L 



to an even stronger approximation. The power lost per unit length will 

be approximated as the smaller of P . . and P , „ 

UJ (volume) tu (surface) 

Using Eqs. (5.29) and (5. 30) we find that the plasma is a trans- 
parent volume radiator for all u) such that 2 R s 1 is true, or 
equivalently, such that 


e (T, p) <; 

uu 


21 


(5.31) 


is satisfied. The transparent radiation losses can be immediately in- 
corporated into the fluid dynamic equation as a loss term. The divergence 
of the plasma flux at any energy uj is simply 


V' B = 4 TT e (T, p) . (5. 32) 

If the opposite inequality is valid, i.e. 

B (T) 

~ p ' S e (T, p) , (5.33) 

2 R uu 

the plasma acts as a blaekbody radiator with effective volume radiation 

2 TT B (T)/R 
UJ 

It is not clear whether or not this loss should be included in the fluid 
mechanics as an effective volume loss. Even though it is a loss from 
the material, the loss occurs mainly at the surface; particularly when 
2^ R » 1 is satisfied. Then most of the plasma in our model is un- 
affected since the temperature is assumed to be held fixed throughout 
the radius. In a more realistic approximation, radial temperature 
gradients must exist which then lead to heat transfer from the hotter 
interior to the cooler surface. The heat transfer is accomplished through 
both thermal conduction and radiative transfer. In this manner the surface 
losses do influence the interior. However a realistic radial profile may 
also result in photon trapping especially for regions of the spectrum 
where is very large. Without additional information concerning the 
radial temperature profile it is difficult to assess accurately the role 
of the surface radiation in fluid equations, especially when the equations 
do not admit radial temperature gradients. This ambiguity suggests that 
two predictions of radiative losses should be made. For an upper limit on 
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the flax which leaves the plasma, both the blaekbody and transparent 
radiation are included. For a reasonable lower limit on radiative loss 
only the transparent radiation and any borderline blaekbody radiation 
(say 2 R < 3) are to be included. 


5.6 APPLYING THE MODEL 


Our criterion is illustrated in Figs, 5, 1 and 5. 3 where, in addi- 
tion to the contributions to the emission coefficient, we have plotted the 
function B^ (T)/(2R) for two values of R of interest. In Fig, 5. 1 where 
the classical approximation to the continuum radiation is shown, the Low 
energy region of the spectrum is seen to lie below the curves ©^(TJ/ZR 
for both radii. Therefore the transparent requirement is satisfied. 

The vacuum ultraviolet portion (a) ^ 13. 6 eV) lies slightly below 

(T)/2R for R = 0, 05 mm but well above for R = 1 cm. Hence the 
VUV is barely transparent for the small radius and strongly black for the 
large radius. Four distinct Lyman lines are shown in Fig. 5, 3, The 
three lowest lines are strongly black at their peak but they possess wide 
transparent wings. The higher lines tend to mimic an extension of the 
continuum. In our calculations of power lost we have lumped the higher 
Lyman lines, Lg and above, with the true continuum, insofar as deter- 
mining whether they are black, marginally black, or transparent* To 
determine the power lost in the transparent limit for the region of 
0 } £ 13.6 eV we have used the quantum -mechanical expression {5. 17) 
for the continuum., and the sum of line strengths for n ^ 4 for the high 
lying lines. For the three isolated lines the centers have been treated 
as black. Fortunately the transition between the black center and the 
transparent wing occurs far enough from the line center that the 
asymptotic Holtsmark profile is adequate to describe the line. The 
transition point Ad) where black turns to transparent can be determined 
by finding Au) such that as found from Eq. (5,25) satisfies the equality 
in Eq. (5. 31 ), 


The total power per ste radian radiated by the wings is thus found 
by integration to be 

C 


line wings 


£ 

3 


nn 


(dr 


3/2 


,/ , 


e duj 
uu 


(5. 34) 


where is the reduced wavelength corresponding to the transition from 
black to transparent. 


The only other line which is strongly reabsorbed in our calcula- 
tions is the first B aimer line The self absorption is important 
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only at high pressures and large radii. The asymptotic formulas are 
not sufficient to estimate the reabsorption, thus the detailed profiles 
of Griem26 must be used. 

The total power emitted per unit volume, including all transparent 
radiation and an effective vacuum ultraviolet radiation whenever it is 
marginally black, is illustrated inFigs. 5.5 and 5. 6 as a function of 
temperature for two values of pressure and radius.* In the numerical 
fluid mechanical calculations of Section 7 , the radiation loss terms are 
adequately approximated as constant factors multiplying one of the semi- 
classical formulas (5,14) or (5.15). These approximations are also in- 
dicated in Figs. 5. 5 and 5.6, All the formulas are of the following 
forms 

-® /T 

p t ,8.6*i<T 9 c, (jt- + c 2 ) x. „ AC e I , JL_) < 5 ., S ) 

I I cm 

-3 

where is in cm . The constants Cj and C 2 used for various 

pressures and radii are converted to a different form, in mks units, in 
Eq. (7.38), and then listed in Table 7.3. 

The calculations listed above are reasonable estimates of the 
effective Loss term per unit volume from the bulk of the hot plasma. 
Radiation which is subject to strong reabsorption has not been included 
in these estimates. The blackbody surface contribution can be easily 
included by Integrating tt (T) over the region of the spectrum which 
satisfies our criterion for blackness. In order to compare the magnitude 
of the black and transparent losses, either the blackbody losses must be 
rewritten as effective volume losses, or the transparent losses must he 
transformed into an effective surface flux. In Table 5.4 we have com- 
pared the black and transparent effective surface fluxes for two situations 
of interest. 

The transparent radiation accounts for the major portion of the 
total power lost. However the blackbody contribution is still a significant 
fraction of the total loss. Furthermore the fraction of the total emission 
which is blackbody radiation generally increases as the radius becomes 


Radiation which is black is not included because of uncertainty in the 
radial temperature gradients, as discussed below Eq. (5.33). 
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Transparent power lost including lines 
2 P,J C ( E q 5.15) numerical fit 





Transparent power lost, including lines 
^ Pq C (Eq5.l5) (numerical fit) 


Hydrogen 

R =.05 cm 
p= 10 atm 


T ( 10 K ) 

Power Radiated from Hydrogen (10 atm, R — .0 5 em) Showing 
Numerical Results and the Analytic Approximation 




TABLE 5.4 

EFFECTIVE SURFACE RADIATION LOSSES 


LASER 10 kW 5 MW 


Pressure 

3 atm 

3 atm 

Temperature 

14000 K 

14000 K 

Radius 

.056 cm 

1,25 cm 

Transparent 

2 

1579 watts /cm 

6507 watts /cm 

Black 

2 

462 watts /cm 

33 56 watts /cm 

Total 

2 

2041 watts /cm 

9863 watts /cm' 


larger. For accurate results at Larger radii it is necessary to include 
the blackbody radiation properly, but this is difficult unless one has 
knowledge of the radial temperature profile, as discussed below Eq. (5, 33). 

5.7 SUMMARY 

Having analyzed the gross behavior of the spectral emission co- 
efficient as a function of w, T and p, we have obtained estimates of the 
volume power loss due to transparent radiation, and of the black body 
surface radiation. These estimates are based upon a simple model - 
one in which the radiation lost per unit volume from the core at given T 
and p is assumed to equal the radiation lost per unit volume from an 
infinite cylinder of equivalent T, p and R. The reabsorption effects 
included are those of the cylinder at uniform temperature and pressure. 
Only two extremes of reabsorpfeion effects are recognized in the final 
calculations of power lost: the radiation of a given frequency is treated 

as being either transparent or black. 

Our calculations show that radiation losses rise rapidly with 
temperature, peaking at approximately 18,00'G K where ionization is 
complete. It is important therefore to know the operating temperature 
accurately, A small error in determining the operating temperature can 
lead to a large increase in radiation lost, Since the operating temperature 
is intimately related to the propagation mechanism, it is also important 
to understand thoroughly how LSC waves propagate. 

Black body surface radiation is less important than transparent 
radiative losses for the configurations studied. However, the black body 
radiation is not negligible, and its importance increases with plasma 
diameter. Future work should include modeling of the radial temperature 
gradients so that better estimates of the surface radiation logj: can be 
made. 


- 49 - 


6. ESTIMATE OF LASER SUPPORTED COMBUSTION WAVE PROPERTIES 


To enable the flowing hydrogen core of the laser -heated rocket to 
absorb laser energy, it mast be heated to a temperature where thermal 
ionization takes place, i.e., to temperatures of order 10, 000 K* In steady 
state the mechanism for providing this heating is a laser supported combustion 
(LSC) wave, which is a zone of rapid temperature increase from ambient 
temperature to temperatures where absorption of the Laser energy can take 
place. Laser energy flows through this zone from the cold side and is 
absorbed on the hot side* The hot plasma transmits energy forward to the 
cold gas by conduction and by radiation emission from the plasma, this 
providing the energy necessary to heat the cold gas up to temperatures where 
it can absorb laser energy* Thus the wave is supported by the laser energy, 
which indirectly converts cold flowing gas to hot flowing gas, using laser 
absorption in the hot gas, plasma emission, and conduction and plasma 
radiation absorption in the cold gas as mechanisms. 

Such waves have been observed in air, in a non -steady geometry, 
by focusing laser intensities of 10§ to 10? W/cm^ on solid targets. 

Evaporation of the targets provides electrons which initiate absorption of 
laser energy in the gas in front of the target, and a temperature wave is 
observed to move away from the target toward the laser beam., at subsonic 
speeds. This wave is an LSC wave moving relative to the still air. It is 
observed to move at a definite speed, which depends on the laser intensity* 

{At intensities greater than 10^ W/cm^, the rapid energy deposition in- 
creases the pressure greatly and caused shock waves to precede the heated 
zone, a phenomenon called a Laser -supported detonation (LSD) wave ). 

The laser-heated rocket flow is conceived to begin with a 
stationary LSC wave providing the mechanism to heat the flowing hydrogen. 

The wave is made stationary by flowing parallel to the laser beam at just 
the speed with which the wave would want to propagate towards the beam. 

This wave forms the initial condition for the flow in the plasma core, so 
it is necessary to know its speed for a given laser intensity, to establish 
the rate of mass flow into the core. 

Experimentally, LSC waves in air have been observed by the non- 
steady method described above and curves of intensity vs. speed derived 
(Ref. 28). The 0 re|icaL^y, a number of attempts have been made to model 
air LSC waves ^ , although their success in comparing with experi- 

ment has been limited. There is no information, either experimental or 
theoretical, on LSC waves in hydrogen. Since radiation from the gas plays 


-SI- 

FAQ® BLAHS NOT FSMU ' 


a crucial role in LSC wave properties, and since radiation properties of 
hydrogen are quite different from those of air, one cannot easily extrapolate 
from air to hydrogen. The lack of reliable information on the intensity vs, 
speed relation of bSC waves in hydrogen must be considered one of the 
major uncertainties in the study of the flowing core laser-heated rocket. 


Nevertheless, it was necessary to estimate the intensity vs, 
speed curve for hydrogen in order to provide initial conditions for the 
plasma core streamtube calculation* For this purpose, the simple model 
of Jackson and Nielsen^ was adapted to hydrogen* In this mod'el we needed 
to characterize the high temperature emission and low temperature absorp- 
tion of hydrogen* For the latter, a classical continuum absorption coefficient 
was assumed. For the former, the radiated power was taken to be the 
continuum radiated power per unit volume, as found by Yos in Ref. 32. 

With these assumptions the Jackson-Nielsen model was used to obtain the 
intensity vs. speed curve for hydrogen at 1 atm and 3 atm. ^ 


The results are shown as the solid lines on Fig. 6, 1. Also shown 
there is the Jackson- Nielsen result for air at 1 atm as the dashed line, and 
the data for 1 atm air of Ref. 2 8, It is seen that the hydrogen estimates do 
not differ much from the air theory in the range 4E5 to 1E6 W/ em^ ( Also 
evident is the disagreement between the air theory and experiment, with 
the theory about a factor of 3 lower than the experimental results. 


In spite of the uncertainty, a point near the 3 atm line at 3. 67E5 
W/cm^ (3. 67E9 W/m^) and 1340 em/s (13.4 m/s) was used for the initial 
conditions of the core calculation. This corresponds for T Q “ 300 K to a 
mass flow per unit area of 0.3316 g/cin^s (3.316 kg/m -s). 

For the 10 atm core calculations, it was assumed that the inten- 
sity vs, mass flow per unit area was the same as for 3 atm; this reduces 
the speed to 400 cm/s (4* 0 m/s). 


These estimates were made before the work described in See, 5 was 
undertaken to define hydrogen radiation emission. That work could be 
used to develop a much more sophisticated model of L.$C waves in hydrogen, 
which avoids most of the assumptions inherent in the model of Jackson and 
Nielsen, and which uses a much more accurate description of radiation 
emission and absorption than is possible with the formulas of Ref, 32. 
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7. FLOWING CORE STREAMTUBE MODEL 


7. 1 INTRODUCTION 

The central streamtube of the flowing-core laser-heated rocket 
combines a number of physical phenomena. It absorbs laser radiation, 
it dissociates and ionizes, it emits plasma radiation, and it flows through 
a variable area channel. All these phenomena must be included in a model 
to calculate the properties of the streamtube. The analysis described in 
the final report on the previous contract (Ref. 1) included laser radiation 
absorption in a variable area channel, but treated the gas as perfect, with 
constant specific heats, and did not include loss of energy by plasma radia- 
tion. Hydrogen is fully dissociated at about 800'OK (with dissociation energy 
2. lb x 10 J/kg) and is 50% ionized between 16,000 and 17.000K (with 
ionization energy 13. 1 x 10® J/kg), while the total enthalpy change from 
300K to 8, 000K is 3. 76 x 10 J/kg and from 8, 000 to 16, O'Q0K is 8, 37 x 10^ 
J/kg. Thus it is clear that the dissociation and ionization processes are 
important energy sinks, and hydrogen should be treated as a real gas. The 
loss of energy by plasma radiation emission both reduces the energy in the 
streamtube and defines an energy loss against which the walls must be 
protected, so it is also important to include this physical phenomena. 

A model of the flowing-core streamtube was constructed which 
included variable area, one-dimensional flow of dissociating and ionizing 
hydrogen in chemical equilibrium, with laser energy absorption by inverse 
Bremsstrahlung, and plasma radiation emission. This model was reduced 
to the solution of three first order differential equations which were solved 
by straightforward integration on a computer. This section will describe 
the model, and give the results of some calculations for configurations of 
interest. 

7. 2 BASIC EQUATIONS 

The one -dimensional mass conservation equation for a gas of 
density 0, speed u in a channel of area A gives the mass flow rate as 


m = p uA, 


A = TT R 2 


(7. 1) 
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The momentum conservation equation relates the pressure gradient to the 
velocity gradient by 



+ 0 u 



(7.2) 


The energy equation has laser absorption terms and plasma emission terms. 
The conservation of energy requires the divergence of the total energy flux 
to vanish. This total energy flux has three parts in the present problem. 

One is the total enthalpy flux of the gas, p uAH, where H is the total enthalpy. 
The second is the laser energy flux, AI, where I is the laser intensity per 
unit area, and the third is the plasma emission radiant energy flux vector 
jh The energy equation is then 


d 

dx 


(p u AH + 



2rr v V 


S dr = 0 


(7.3) 


The laser energy is absorbed with absorption coefficient accord- 
ing to 


d At 
dx 


k AI 

Jj 


(7.4) 


who 


se solution is, with { denoting conditions at x = 0, 


x 

AI = (AI)j e T , T = J* dx 

0 


(7.5) 


Use of Eqs. (7. 4), (7, 5) and {7. 1) in (7. 3) allows the energy equation to be 
written as 


dff __ <|h udu _ dQ 

dx dx dx dx 




k T (AI). e" T 


R 

J 2wr V ' 


A 


I dr 


(7.6) 



This is written so the laser absorption term is recognized as a heat addition, 
and the plasma emission term as a heat loss. 

The problem is then one of heat addition to a one-dimensional nozzle 
flow, with the heat addition term coupled to the flow, since and V ♦ J> 
depend on the local properties of the gas. 

The conventional way to solve this problem, as described in Ref. 1, 
is to take the area A(x) as given, and solve for dp, dia and dp from {7. 1), 

(7. 2) and (7. 6), The results show the sonic singularity at M - 1, which 
makes for difficulties in the numerical integration of the equations when 
there is heat addition or real gas effects, both of which are present here. 

An alternative method of solution was therefore used, which led to simple 
numerical solutions. 

Ii ead of taking the area as given, we take the velocity distribution 
u(x) as known. This leaves the thermodynamic variables p, p, T to be 
found from (7. 2) and (7. 6) and the equation of state. Then A(.x) is found 
from (7. 1), since the mass flow rate m is given. This is a perfectly 
acceptable method of solving the problem, as long as one chooses u(x) so 
the resulting area distributions are realistic, and this can be accomplished 
with a little experience. 


Equations (7, 2) and (7. 6) can be written easily in terms of dp and 
dT by using the equation for the enthalpy of equilibrium hydrogen in the 
form 


h = h (p, T) 



(7.8) 


If this is inserted into (7.6), and {7. 2) used to eliminate dp, (7, 2) 
and (7. 6) become 


When the thermodynamic properties of hydrogen and the radiation properties 
which are contained in dQ are defined, the specification of u(x) enables a 
straightforward integration of these two equations from any given initial 
condition. From their solution and the gas law, the radius of the stream- 
tube is obtained from (7, l) as 

R = (m/TT o u) 1/2 (7. 10) 

Actually the emission term in (7. 7) involves R, so (7, 10) is used for R 
during the course of the integration, 

7, 3 THERMODYNAMIC PROPERTIES OF EQUILIBRIUM HYDROGEN 


We need the derivatives of enthalpy with respect to p and T, the 
equation of state which relates q and T, and the composition of the mixture 
of hydrogen molecules, atoms, ions and electrons. The pressure deriva- 
tive of h can be expressed in terms of the equation of state by a classical 
use of reciprocity relations between the derivatives of h and entropy, which 
are in turn obtained from the second law. These lead to 


so that 



(7,11) 


We characterize the hydrogen as a mixture of perfect gases whose 
components have number densities: molecules n^, atoms n^, ions nj and 
electrons mg. We do not expect to get to a temperature where other species 
are significant, nor are there multiply ionized atoms, so nj = mg. As 
reaction coordinates we use the fraction of molecules dissociated P and the 
fraction of atoms ionized, a : 


* n A + n i 

? ’ n M ' (:! A ' "r 1 '- ' a " "A + *1 

In terms of the original number of particles per unit volume. 


(7. 12) 


m - n., + ( n A + n ) / 2 
o M A I' 


(7. 13) 


the number density of the species are 


n = n Q (1 - 0) . n A = 2n 0 S(l - a) , ^ = n^, = 2 h q @ a (7. 14) 
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and the ear responding partial pressure are found from 


p. = n.kT (7.15) 

where k is the Boltzmann constant. The partial mass densities are, in 
terms of the particle masses rm, 

P M = n M m M ’ P A "A ®A ' °I = "I "V °E *1 “ E < 7 ' 16 > 
and the total number and mass densities are 

n = n Q [1 + fi (1 2a)] . p = n Q ( 7 . 17 ) 

Whenever convenient we will use the facts that to an excellent approximation 
for hydrogen, rn^ = rrij = m^/2. 

The equation of state is the sum of the partial pressures. 


p = D Z R M T , Z = 1 + 8 (1 + 2a) , R M = k/m M (7.18) 

where Z is the compressibility factor and the gas constant for the mole- 
cules. We have now expressed p in terms of p, f through a and I by (7. 18). 
Since p and T are the dependent variables in the basic differential equations 
(7. ?), we need q, 8 in terms of T, p. 

Since we deal with hydrogen in thermochemical equilibrium, a and 
8 are determined by the Law of Mass Action in terms of the thermodynamic 
variables ©f the mixture. There are two reactions occuring, dissociation 
and ionization: 

H 2 ^ 2H H — H + + E‘ 

For each the equilibrium constant provides a relation between the 
partial pressures and the partition functions of the species, involving the 
heat of reaction. These relations are 
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where h is Planck's constant, Q v and Q are the vibrational and rotational 
partitioi? functions of K , are the electronic partition functions of the 

species, h^ 0 is the heat of dissociation of a molecule per unit mass and 
hj° is the heat of ionization of an atom per unit mass. From (7. 14), (7. 15), 
(7. 17) and {7. 18) the partial pressure ratios are 


£a 4B 2 H - l)' 1 ' . 48 z (1 - a) 2 J>. IV. 20a) 

p M - t -f % kT ‘ 1-6 2 
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(7.20b) 


The rotational partition function Q may be approximated by T/2 0 where 
® r = 87. 62K for H 2 and the factor of 2 comes from the symmetry of the 
molecule. This is a valid approximation for T » © j it is about 10% low 
at 300K but only 0. 5% low at 1Q00K. The vibrational partition function is 

Q v = ^1 - e j where © v = 5983K for H^. The electronic partition 

function for the molecule may be replaced by its ground- state statistical 
weight 1, since its first excited state lies very high (132, O0QK). For the 
electron Q e jj* * 2 because of the two electron spin states, and for the ion, 
there are no internal states, Q y. = 1. 




For the atom 


Q elA =2S 1 eXf> { " h i° m I /kT ) 


(7.21) 


17 

V' - 2 o 2 

7 , J exp (hj m^/kT j ) 


The first term is the usual ground- state weight 2. The subsequent terms 
are needed here, even though hj m . /k - 158, O0OK, because we may be 
interested in hydrogen atoms up to 20, @@0K. The upper limit cut-off of 
the sum can be obtained in various ways, and may depend on temperature. 
Here we took it constant at 17, since the effect of the sum is rather small. 
If {7. 20) and (7, 21) are used in {7. 19) we have 
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These two relations express a, P in terms of p, T, as desired, and could be 
used* However, they are algebraically complicated, and difficult to solve 
for a, P. For hydrogen they can be considerable simplified with little loss 
in accuracy. Hydrogen in equilibrium dissociates almost completely before 
it ionizes significantly* We can thus separate the two reactions and take 
a - 0 in the first one while I goes from 0 to 1. Then in the second, we 
can take P = 1 while a increases from zero. These two are separated by 
a temperature T below which there are only molecules and atoms, and 
above which there are only atoms, ions and electrons. Additionally, in the 
non^ ionizing region, the temperature is too low for the higher electronic 
states of atom to contribute, so the squared term on the right of (7* 22a) 
can be taken a$ unity. Then we have 
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' 2TT m E k \ 3/2 kT 5/2 

V T > MT 7 

v h p 
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These are both linear in (3 ora , and yield 


T < T* : a = 0, Z = 1 + 0 = (1 + p/fjj)" 1 ^ 2 (7.24a) 

o, 1/2 

T>T^: B = 1, Z = 2 (1 + a), a = (1 + p/fj) ' (7.24b) 

which are very simple expressions for gl, $ in terms of p and T, the pri- 
mary dependent variables. 

From (11), we need the derivative of In p with In T at constant 
pressure, which from the gas law (7. 18) is 


/ S In p \ , / I In z\ 

ln T/ p = ' ‘W*T/ p 


(7.25) 


The derivate of Z is easily found from (7. 24) in terms of the derivatives 
of a, 8, which are obtained by differentiating (7. 24) and using the definitions 
of f^, f^ given in (7. 2 3). The re suits are 
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T < T : (d a/a T) = 0, i|r = 

p y v 


= (v t )/C vt - 




I **■ ) = MI- 8 2 ) + !fi. . * 1 

In T J 2 \JL T vj 


(7.26a) 
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(7. 26b) 


The remaining thermodynamic property needed is Cp, the tempera- 
ture derivative of the enthalpy. The translation specific heats of the atoms, 
ions and electrons are 5 k/2 per particle. This is the total contribution for 
the ions and electrons, which thus have a specific heat per unit mass of 


v * ,k/2ra i • c ,E*’ k/,m E* V ”l /m B 


(7.27) 


Their enthapies are the integrals with respect to T, plus the dissociation 
energy h^ 0 and the ionization energy hi° for the ions 


h = c _ T + h° + h_° , h = c 
X pi D I E 


T 


(7. 28) 
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The atoms have a translational contribution 5kT/2m^, carry the dissocia- 
tion energy hj-j°, and also have a contribution from the electronic states 
derived from Q e ]^ of (7.21) by 


h -ill! 3 10 °eIA (j 3 

elA m A IT I \ Sj 


where has been defined in (7.26b). The resulting atom enthalpy is 


h A = + h n 0 + h T° t 1 -T 

A 2 m^ D I \ Sj 


(7. 29) 


Finally the molecules have a translational contribution 5kT/2m^, a 
rotational contribution 2 k T/2 for the two rotational degrees of freedom, 
corresponding to Q = T/2 ® r , and a vibrational contribution found from 

-© /T 

_ , v 

Q = 1 - e as 


h = 


2 a In Q . 
k T v _ k T 


v m 5 T 

M 




where tji is defined in (7.26a). Thus the molecular enthalpy is 


h . 2U ♦ ii , 

m 2 cn. , m. v 
M M 


(7.30) 


wei 


The enthalpy of the mixture is the sum of the species enthalpies 
i by their mass fractions in the mixtures 


h = 52 p /@ 


(7.31) 


The mass fractions are found in terms of a and [3 by (7. 14), (7. 16) and 
(7.17): 

P m °a P I °E m E 

_JE = 1 _ b . -A. = e (1 _ a), -=Ba, T =Ba — (7,32) 
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Using (7.28), (7.29), (7.30) and (7. 32) in the enthalpy expression (7. 31), 
we find the enthalpy for this mixture as 


h = c pI T | (B (1 +a) + (1 - 8) (0. 7 + 0.2 * v ) 

+ 0.4 8 [© d / 2T + (a+ (1 - a) t ftlA ) Ij/T J j 


(7.33) 


c ... = 5 k/2 m , * , = 1 - (S,/S_) 

pi I elA 1 2 


The equilibrium is obtained by differentiating h with respect to 
T at constant p: 
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(1 - 0) (0* 7> 0. 2 cp ) + 8 (1 + a) + 0.4 g (1 - a) 
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+■ Gt + 3,2 ^ + 0*4 



elAJ 



(7* 34) 


cp = (8 /T) 2 

V V 


yV( c 'v /T 
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17 

S 4 = ^ j" 2 ex P C®j/T j 2 ) 

3 = 1 


Equations (7.11), (7. 18), (7.23), (7.24), (7.25), (7. 26) and (7.34) 
express ©, (bh/Sp)^ and e in terms of p, T, and thus express all the 
unknowns except dQ in the oa sic differential equations (7.9) in a form 
usable for numerical integration 
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A cheek on the thermodynamic properties used can be obtained by 
comparing some calculations made using the present formulas with those 
of Patch {Ref. 33). (The enthalpies of Ref. 33 have had h D ° = 2. 16 E 5 j/g 
added, since in that reference the H atom is taken as the reference for 
enthalpy, while here we have taken the H? molecule. ) The comparison is 
made in Table 7. 1 for 8, a and h, with T - 8000K. Comparison shows the 
errors to be mostly a few percent, compared to the very elaborate scheme 
of Ref. 33. The worst errors are at high values of T and p. The applica- 
tions in the present report do not exceed 16,00OK at 10 atm, and are all 
above 8, 000K, so the approximations appear satisfactory. The major 
source of error is probably the use of 17 terms in the sums of Sj and 
instead of a variable cut-off. 

7.4 ABSORPTION OF LASER ENERGY 

The first term in the heat addition expression, Eq. (7. 7), gives the 
absorption of the incident laser energy by the plasma. It is dependent on 
k^» the absorption coefficient of the plasma. Based on the studies of Ref. I, 
we use as the absorption mechanism inverse Br emsstrahlung, both electron- 
ion and electron-neutral. The expressions for k^ have been given in Ref. 1. 
The electron-ion absorption is given in Eq. (2. 1 ) of Ref. 1 as 

k LEI = S EI "e ” I [»*P (lyv/kT) - *] 


where j- s the electron-ion absorption cross-section and v is the laser 
frequency (for 10.6pm in the present case). The cross-section is given 
in Eq. (2.2) of Ref. las 


a 


El 


4 

3 



1. 63 x 10 ^ , 5, 

... 1/2 - (cm > 

T 


where Ze is the ion charge, e is the electron charge and e is the speed of 
light. Combining these two equations we find 


1.63 x 10' 32 , 13S7/T ,, , -I, 

k EEi ' "T7F "e ‘i <e - *> (cm > 


where T is in K and n , n_ are in criT 


TABLE 7. 1 


COMPARISON OF tt. g, h WITH REFERENCE 33 


r (atm) 

T (K) 

_L 

B (Ref. 2) 

% Error 

1 

2000 

8. 34E-4 

8. 10E-4 

3. 0 

10 

2000 

2.64E-4 

2. 56E-4 

3. 1 

1 

4000 

0, 645 

0. 622 

3. 8 

10 

4000 

' 0.258 

0.243 

6. 0 

1 

6000 

0. 99389 

0. 99251 

0. 14 

10 

6000 

0. 9435 

0. 9319 

1. 2 

1 

8000 

0. 99994 

0.99926 

7.3 

10 

8000 

0. 9945 

0. 9927 

1.9 

p (atm) 

T CK) 

a . 

0. (Ref. 2) 

Error 

1 

10, 000 

0. 2126E-1 

0.2192E-1 

-3. 0 

10 

10, 000 

0. 672E-4 

0. 706E-4 

-4. 8 

1 

14, 000 

0. 2922 

0. 3050 

-4.2 . 

3 

14, 000 

0. 174 

0. 185 

-5.9 

10 

14, 000 

0. 0962 

0. 105 

-8.4 

1 

20, 000 

0. 895 

0. 922 

-2. 9 

10 

20,000 

0. 535 

0.651 

-17. 8 

p (atm) 

X .(E) 

k 11/ R) 

h (j/g) (Ref, 2) 

% Error 

3 

12,000 

5.604E5 

5. 57E5 

0. 5 

10 

12,000 

5. 199E5 

5. 16E5 

1. 1 

3 

14,000 

8. 123E5 

8. 03E5 

0. 8 

10 

14,000 

6. 908E5 

6. 74E5 

2. 5 


The electron and ion number densities are related to the density by 
(7. 14) and (7. 17) as 


n 


n T - 2 (So. 0 /m.. 
I M 


The above expression for k T _ -ir . T is a semi- classical result. To 
account for quantum - mechanical effects we multiply by a Gaunt factor of 
1. 5, which is an average value for the range of interest. The final expres- 
sion for the electron- ion inverse Bremsstrahlung absorption coefficient 
is then 


k = 8. 7 x 10 
LEI 


2 2 2 
a | p 


. 1357/ T ,, 
(e - 1) 


{7. 35a) 


a - 15 for egs units, a = 11 for mks units. 

Of course, this expression is only valid where there are ions, a =£ 0, so 
in our approximation for the chemistry, it is to be used only above T * , 
where 1=1. 

The absorption coefficient for electron-neutral inverse Bremsstrah- 
lung is also given in Ref. 1, Eq. (2. 6), as 

k LEN = °EN "E "N <-V' kT >] 


where Qjgjq is the absorption cross- section and njq is the neutral number 
density. Equation (2. 8) of Ref. 1 gives 


Q 


EN 


2, 96 x 10" 45 T 
I - exp (-h v/kT) 



where £ may be approximated by the first term of Eq. (2. 9) of Ref. 1, 
when VT/ 158, 000 is small, as 

1/2 

?= 4. 862 ('kT/m h®) 
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Since for hydrogen nij hj°/k = 158, 000 K, then we have 


= 1. 6 x 10 38 JW n n exp (-0. 0122 Jt ) (cm" 1 ) 


- 1 . 


LEM 


E N 


-3 


where ng, n^ are in cm - . 


The number densities are found from (14) and (17), where we 
recognize that our chemical model has 8 = 1 whenever there are electrons 
present, so 

“s** #e/ "V n N = n A = 2 11 " *> Dm M • 

Then the final expression for kggg is 


'LEN 


= 5. 7 x 


io b a 


(1 - a) p 


2 


exp (-0, 0122 


(7. 35b) 


b ~ 9 for cgs units b = 5 for mks units 

The expression for kg to be used in the first term of (7. 7), and in 
the definition of optical depth T, Eq. (7. 5), is the sum 


K L k LEI + k LEN 


(7,36) 


A plot of kjjgj (dashed lines) and (solid lines) is given in Fig. 7. 1 as 

kg/p vs T. It can be seen that electron-ion absorption dominates at temp- 
eratures above 11, 00 OK at 10 atm, and above 9, 000K at 1 atm. The sum 
kg follows the kggj curves closely above 12, 000K. 

7. 5 PLASMA RADIATION EMISSION 


The second term in dQ, Eq. (7. 7), is the loss of energy caused by 
emission of radiation from the plasma 
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2tt r V • 


1 dr 
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The model used for this term is discussed in the section of this report 
entitled Radiative Emission. The model is one of a uniform cylinder of 
plasma, so the emission at any station is expressed in terms of a volume 
emission and the cross-sectional area 

R 

y* 2TT r V ’ J dr = TT R 2 F t ^ 

0 


The general form of as deduced in the section on Radiative 

Emission, from Eq. (5.32), is 


P T = 8.6 x 10" 9 C x 
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AG 


_X 

0 


~® r /T 



- 3 

where n^Q is the atom ground state population^ in cm , and the constants 
and depend on the pressure level and the relative importance of line 
and continuum radiation. The relation between n^^ and n ^ is through the 
electronic partition function Q e ^ of (7-21), whose first term relates to the 
ground state. 


n 


n ^ = n 
AG A 


'elA Sj exp {- h^° m^/kT) 


When n^ is expressed in terms of a> £ and p from (7, 14) and (7, 17), the 
expression for is, in cnks units. 


p - * i 

■ T ^ 1 


'•'*sg • =,) f 


T e.(l - a) Q 


s 


1 


6 ) 


(7. 38) 


The values of the constants C^, depend on the pressure level and laser 
power, as described in the section on Radiative Emission from Hydrogen, 
The particular values used for the cases calculated will be given when the 
results of the calculations are discussed. 


7. 6 SPECIFICATION OF VELOCITY DISTRIBUTION 


As described above, in the subsection on Basie Equations, the 
velocity distribution in the plasma core str eamtufee is chosen as the indepen- 
dent parameter rather than the area distribution, since it greatly simplifies 
the calculation of the flow, by avoiding the explicit appearance of the singu- 
larity at the station at which M = 1. 
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Several forms of u were tried, to see if reasonable nozzle shales were 
attained. The simplest and most satisfactory was a logarithmic distribution. 

The calculations were done with the optical depth t, rather than 
distance x, as the independent variable. This choice eliminates the explicit 
appearance of from the differential equations when the plasma emission 
term e is not taken into account, as can easily be seen from Eqs. (7.5), (7.7) 
and 7.9). It was used in Ref. 1 , where e was not considered. In the present 
case where e was included, there is no simplification afforded by using T, 
since it will require the explicit appearance of k^ in the plasma emission 
term, but it was used to provide some continuity with Ref. 1. However, in 
the future, it probably would be better to use x as the independent variable. 

In terms of T, the u distribution was taken as 



which goes from Uj at T = 0 (x = 0), the initial station, to u f at the final 
station T £, This is a linear variation of In u with T. The nozzle shapes 
obtained with this can be quite reasonable, as will be seen when the results 
of the calculations are presented. 

The derivative is easily obtained from 


as 


In u = In + (T / T^) In (m^/u^) 


du 

dT 


u In (u f /u ) 
T 


(7. 39b) 


7. 7 SUMMARY OF THE MODEL 


All the quantities necessary for the integration of the differential 
equations (7.9) have now been defined. In terms of T as the independent 
variable, (7. 9) becomes, using (7.5), (7.7), (7.10), (7.11), (7.18), (7.25) 
and (7. 37): 



" £> u 
R m ZT 


du 

dt 


(7. 40a) 
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(7. 40b) 


dx 

dT 


R = 


m R. ZT 
M 

pm 


1/2 


(7. 40e) 


Z is given in terms of a, [3 by (7. 24), and those equations also give g» $ 
in terms of p, T. Equation (7.34) gives c^, while the derivatives of a, I 
and Z are in (7.26). With m and du/dT from (7, 39), a forward integration 
of the three differential equations can be performed from a starting station 
T = 0 at which x = 0, T = T\, p = p^, u = u^ to a final station T^, where 
u = Uf. The mass flow rate constant m is determined by the given initial 
conditions and the initial radius R^ as 


m = TT R. p. u. = A. p. u. 

ill ill 


(7=41) 


which shows that only R/Rj is determined by the solution. Of course, if 
the laser power P is given, R^ is related to the incident laser intensity by 


I = P/TT R. 
© l 


(7.42) 


The incident laser intensity I is not the same as the initial intensity Ij, 
because energy has been added to the incident gas in the LSC wave before 
the stream tube calculation is started, and we must be careful to conserve 
energy. If the calculation is begun at T^, p^, u-, the initial total enthalpy 
flux is m (h* 4- u^/2), where h^ is obtained from (7.33). If the gas i^s 
supplied at p 0 , T Q , h e , u q , then its total enthalpy flux is m (h Q 4- u Q / 2), 
and the amount of energy added up to the initial state must have come from 
I„. Thus 


I. = I - m (h. + u. /2 - h - u ^/2) 
i o Ail o o 


(7.43) 
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where is the mass flew per unit area, 


= fh/ A . “ 0. u. = p u . {7 0 44) 

A 1 1 l o o 

If the initial conditions are at a temperature of order 10, O0QK, the reduction 
in 1^, caused by the heating of the incident gas by the LSC wave, is sub- 
stantial. 


The initial conditions for the streamtube are determined by the 
properties of the LSC wave which pre-heats the gas. As is customary with 
LSC wave properties, we assume it causes no pressure change, but only 
temperature and density changes, so p^ - p^. In addition the LSC wave 
propagates against the flow at a speed dependent on laser intensity and pres- 
sure level, so I and p determine m A . If T is also given then o and 
u q = m^/o o u., “ rn^/p-. Finall Y* (7.43) gives 1^, and we have all the 
necessary initial conditions* This scheme shows that the properties of 
the LSC wave, plus a choice of l Qt p Q and T Q are all the conditions neces- 
sary to integrate the streamtube equations, once u^ and are chosen# The 
integration then determines the distribution of all gas properties along the 
streamtube including the plasma emission e, and also determines the 
streamtube shape R/R,. The scale of the streamtube can only be found 
when the laser power is also given, whence R^ follows from (7,42). 
However, it should be noted that some idea of the size of the streamtube, 
is necessary to perform the integration, because its size influences the 
degree of opacity of the plasma, and thus influences the choice of constants 
in the expression for G. 

Although it is not used in the integration, the Mach number is an 
interesting, parameter which can be calculated when p and T are known, 
since they determine the speed of sound a. For a gas in thermochemical 
equilibrium, the second law of thermodynamics and the expression of h as 
a function of p and p shows that 

(bh/^o) 

- - R 

s 0 1 - (bh/dp) 

P 

Conversion of these derivatives to h as a function of T and p, and use of 
(7. 11) shows that 
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The second derivative has already been expressed in terms of Z by (7.25). 
From the gas law (7. 18} it is easy to see that 


dp \ 


R M ZT 


§ In Z 


d In p 
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Thus the Mach number squared is 


, ,2 «* _ » 2 

' a J " T I R M Z 


S In z 
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a In z 
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(7.45) 


in {7. 24) with respect to p gives 


T < T 


* / a In z \ 

: \d In p A 


X > T 


: /a In z \ 

: \d In P L 


(7.26). 
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differentiation of Z, 
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' d| \ 

Ml - B 2 } 

(7. 46a] 

1 "+ B ' 
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1 i 

( dq \ 

a (1 - a 2 ) 
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1 + a 1 

id In p / _ 
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Thus the Mach number can be calculated as the integration progresses. 
7. 8 RESULTS OF CALCULATIONS 


A computer program has been written to integrate Eqs. (7.40), and 
has been used for six cases of varying laser power and incident gas pres- 
sure level. All the eases were done for the same incident laser intensity 
of I = 3. 67 x 10^ w/m^ (3. 67 x 10^ w/cm^) and the same mass flow rate 

o • 7 9 

per unit area, = 3.316 kg/m (0.3316 g/cm ), based on an estimate of 
the properties of LSC waves in hydrogen, as described in the section on 
LSC Wave Properties, In all eases the incident temperature was T = 30QK 


and the initial temperature was T. = 12, OOOK. (It is believed that the effect 
of lowering T\ by a few thousand degrees will not change the major results 
of the calculations. ) 

The parameters varied were the laser power P and the incident gaB 
pressure level p Q = p^. Therefore the initial radius R^ (and area A^, and 
the total mass flow rate is also varied. The calculation matrix is given in 
Table 7.2. 

The range of power varied from the experimental device of 10 kw to 
a 5000 kw device of operational size, with an intermediate size of 100 kw. 
The pressure of 3 atm corresponds to the experimental device, while 10 
atm was chosen to show the effect of operating pressure. The initial 
densities and velocities are determined by p Q , Tj and , as explained 
above. In all cases, the integration was carried out to a final value of 
optical depth = 4, with a velocity there of u^ = 1.47 E4 m/s. The attempt 
was to perform the calculations to the throat of the streamtube, and in all 
cases the minimum area was in fact reached. The choice of = 4 means 
that the laser energy s|ill unabsorbed at that point is down from the initial 
value by a factor of e* =1.8 E-2, so only 2% of the initial laser energy is 
not absorbed. From there one, there is no longer significant heat addition 
to the gas. There are still losses due to radiation. 

The constants C-j and used in Eq. (7.38) to express the plasma 
emission e were taken from the material described in the section on 
Radiative Emission, and were appropriate to the sizes of plasma core for 
each power level. They are given in Table 7, 3. 

The properties of hydrogen needed were taken from standard 
references, and are: 


k = 1. 305 E-2 3 J/K, m A = m = m„/2 = 1. 873 E„27 kg 

I M 

m = 9. 109 E-31 kg, h_ = 6. 6258 E-34 Js 

xC Jr 


9 r = 87. 62 I<, 


0 = 5983 K 

v 


0 = 52, 000 K, 


Rj = 158, 000 K 


The value used for includes the correction for amharmonicity. 
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TABLE 7.2 
CASES CALCULATED 


Case 

P (k\V) 

P q (atm) 

(cm) 

/V (cm^) 

. Z 2 

m(g/cm -s=lQkg/m -s) 

i 

10 

3 

9.31 E-2 

2.72 E-2 

9. 03 E-3 

2 

5000 

3 

2. 08 

1. 36 El 

4. 51 

3 

100 

3 

2. 94 E-l 

2 . 72 E-l 

9. 02 E-2 

4 

10 

10 

9.31 E-2 

2.72 E-2 

9. 03 E-3 

5 

100 

10 

2.94 E-l 

2.72 E-l 

9. 02 E-2 

6 

5000 

10 

2.08 

1.36 El 

4. 51 




TABLE 

7. 3 




CONSTANTS IN PLASMA EMISSION 


Case 

F (kW) 

p (atm) 
o 

C 3 (W/m 3 ) 

c z 

1 

10 

3 

2 

2.4 

2 

5000 

3 

2 

0.4 

3 

100 

3 

1 

2.4 

4 

10 

10 

4/3 

2.4 

5 

100 

10 

3/4 

2.4 

6 

5000 

10 

2.2 

0.4 
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The integration of the differential equations was performed at steps 
in T of 0. 05, so there were 80 steps from t = 0 to 4, 

The profiles of T and R/R^ vs x are given in Figs, 7. 2 and 7. 3, 
Figures 7.2a and 7.2b are the temperature profiles for 3 and 10 atm respec- 
tively. The location of the throat, and the station at which 90% of the initial 
laser energy is absorbed are marked. They all show an increase to a rather 
flat maximum, followed by a drop to the throat, although the throat is still 
at temperatures of 13, 000 to 15, 000 K* The maxima are quite flat, and 
occur near the 90% absorbed station. The larger the power, the larger is 
the temperature at any station, as might be expected. Also, the large the 
power, the shorter the distance to the throat, though this effect is not great. 
These distances vary from 3.2 to 4* 0 cm for 3 atm and from 0.49 to 0. 54 
cm for 10 atm. This points up the major effect of pressure level* Because 
the absorption coefficient increases with pressure, the laser energy is 
absorbed in shorter distances at the higher pressure, and the converging 
part of the streamtube is shorter. The 10 atmosphere streamtubes are 
nearly one-tenth the length of the 3 atmosphere streamtubes. 


When one combines the initial radii given in Table 7,2 with the 
throat lengths given in Fig. 7.2, the small size of these devices becomes 
evident* The lengths are determined by the pressure level through the 
absorption coefficient, as just mentioned* The radius is determined by the 
power and the incident intensity I . If I of order 10^ W/cm is required to 
support the L.SC wave, then radii of order 0, 1 to 2 cm will follow* To 
enlarge the size, it is crucial to operate at as small values of I as possible 
which means as near as possible to the threshold intensity of the L SC wave* 
This emphasizes the importance of having a reliable theory to estimate 
the threshold of LSC waves in hydrogen. 


Figures 7,2a and 7. 3b show the normalized radius distribution 
R/Rj for 3 and 10 atm respectively. It can be seen that only small varia- 
tions of shape with power occur, and the shapes seem very reasonable. 

Notice the vertical scale is not physical size, though the horizontal scale 
is, so angles on the figure are not necessarily correct. For the 5 MW 
power R| = 2.08 cm, so the streamtube shown on the 3 atm figure. Fig. 7.3a, 
has very nearly the same scales on both axes. The same is approximately 
true of the 100 kW power on the 10 atm figure, Fig, 7. 3b, where = 0.294 
cm. Other streamtube shapes can be produced by varying the assumed 
velocity distribution, and more experience with calculations would show 
what range of shapes can be achieved* 

Notice that at 10 atm the higher power cases lead to fat, stubby 
streamtubes, such as one which is 2, 08 cm in initial radius and only 0. 5 cm 
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long for 5 MW. This may be undesirable geometry from flow criteria, 
and may indicate that high pressure operation at large powers is xiot a 
useful combination* 

It can be concluded that this rather small program can calculate the 
core streamtube quite easily. 

Other results of the calculations are described in the sections or* 
Radiative Heating and on Performance. 
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8. HEAT LOSSES 


8.1 INTRODUCTION 

The flowing plasma core loses heat by emission of plasma radiation, 
and by convenction of heat to whatever surface is in contact with the hot 
plasma. The walls of the rocket must be protected from this heat load* In 
order to define what method of protection should be used, the size of the 
heat load must be known. This section is devoted to calculation of the heat 
load from both emission and convection. 

The emission has already been used in the plasma core streamtube 
calculation, where it was coupled to the plasma core temperature and 
geometry to provide the energy loss term. The explicit calculation of the 
amount of energy lost and the resulting heat load is easily obtained from 
the results of the core streamtube calculation, as described in Subsection 
8 . 2 . 


No convective heat loss was taken into account in the uor-e calculation. 
In order to decide if this loss is of sufficient size to require it also to be 
coupled to the core calculation, an uncoupled calculation was made by using 
the temperature and pressure distributions found from the core model. 

These conditions were used as external conditions in a hydrogen boundary 
layer calculation to estimate the heat loss tc- a cool surface placed in contact 
with the plasma core. If the resulting energy loss is large compared to the 
core energy, it indicates a coupling of the convective loss and the core flow 
i$ needed. If the loss is small compared to the core energy, such a coupled 
model is not necessary* In addition to an impervious wall boundary layer, 
some cases of gas injection into the boundary layer were calculated to see 
if that was a feasible method of reducing the convective heat loss. These 
hydrogen boundary layer calculations are described in Subsection 8. 3, 

Finally, the possibility of protecting the walls from the radiative 
emission by particles in a buffer gas region surrounding the plasma is 
considered in Subsection 8.4. 

8. 2 RADIATIVE HEAT LOSS 

One of the results of the core streamtube calculation is the radiative 
energy emission from the plasma* At any station, the plasma emission 
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term is TTR P^ # as shown in Eq. (7* 37). The streamtube program 
calculates both R and P'j as function of x, and the quantity TT is 

shown in Fig. 8* la for 3 atm and in Fig. 8. lb for 10 atm. Notice that in 
Fig. 8. la, the 5 MW curve is plotted as TT R 2 P^/lQO, while in Fig. 8. lb 
the 100 kW curve is TTR 2 P t /10 and the 5 MW curve is tt r 2 P T /1O0. 

These curves show that the station emission has a maximum some- 
what earlier than the temperature profiles shown in Fig. 7. 2, because R 
decreases steadily. At the throat, the values are below the initial values 
because of the small radius there. 

Ip order to estimate the heating received by the wall or the buffer 
gas by plasma emission we divide the station emission by the circumference 
2 TTR of the steamtube at the same station, to get the distributions of 


q rad = R R 2 P T /2TTR (8.1) 

The results are shown in Fig. 8. 2a for 3 atm and Fig. 8. 2b for 10 atm. 

(The 5 MW curve on Fig. 8* 2b is 1/10 of q ra( ^), This is not exactly the 
radial heat flux at the station, because the radiation is emitted in all 
directions, not only radially. However, the includsion of the geometrical 
effects necessitates a complicated calculation which is nt warranted at 
the present stage of the theory. For long, narrow streamtube s the approxi- 
mation is a good one since most of the solid angle is radial, and the axial 
energy transport out of one station is compensated by axial transport in 
from other stations. For short, stubby streamtube s like the 5 MW, 10 atm 
one, the approximation is not so good, but these are probably not desirable 
geometries to use. 


The total energy lost to the streamtube by plasma emission up to 
the throat can be obtained by integrating the curves of Fig, 8,1 from x - 0 
to the throat station. The average energy flux at the streamtube surface 
is Likewise obtained by integrating the curves of Fig, 8.2 to tt:e throat, 
and dividing by the throat station x^. 
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Table 8. 1 presents these quantities for the six cases calculated, 
as well as the maximum value of the energy flux 


(q pad ) = (rr R 2 P T / 2nR ) (8.3) 

max 

and the fraction of the laser power radiated. The radiative losses are 
small fractions of the power input, indicating that the incident power is ' 

converted efficiently into hot plasma, which is then available for propul- 
sion. The larger engines are the most efficient; even though the absolute 
amount of power they radiate is larger, the fraction of the incident power 
is smaller. The reverse is true for the energy flux; the larger engines 
have the highest fluxes. This is a surface to volume ratio effect at any one 
pressure. The higher fluxes at the higher pressure are a result of the 
shorter length of the streamtube. The magnitude of the fluxes at the 
higher powers indicate that some form of active heat protection will be 
needed, either by absorbing some of this energy in a buffer gas, or by 
some suitable scheme of cooling the walls. 


§. 3 CONVECTIVE HEAT LOSS 


The hydrogen boundary layers grow on the sides of the core 
streamtube starting at the initial station, with the flow external to the 
boundary layer given by the core properties. For the six cases calculated 
by the core streamtube program., the Reynolds number based on conditions 
in the core at the throat station varies from 2800 for the 10 atm, 5 MW case 
to 8500 for the 3 atm., 10 kW case. This indicates that laminar flow will 
prevail, so a laminar compressible boundary layer model is appropriate. 

Since the primary purpose of this calculation is to estimate the 
heating, and to see if a coupled calculation is necessary, it is appropriate 
to use as simple a model as will include the important physical effects. 

We have therefore assumed that a "local similarity 1 ' solution will be 
satisfactory. That is, at any station., the boundary layer is calculated as 
if it was a similar boundary layer with the local external properties as 
parameters. This reduces the numerical problem to one of ordinary, 
rather than two - variable partial, differential equations, which is a great 
simplification. This approximation has proved to be very useful in 
calculating heating from high temperature gases , and can be /ound 
discussed in books *n boundary layer theory. 

Although many eases of similar boundary layers have been calculated, 
none known to us is useful for the present ease of very high temperature 
hydrogen, because they do not include the transport properties of hydrogen 
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TABLE 8. 1 


RADIATION HEAT LOSSES 



p =3 atm 
o 

Power (kW) 

.... . . .... ..... 

10 100 5000 

Throat Station x (cm) 
tn 

3. 99 3.46 3.24 

Power Loss Q , (kW) 
raa 

2.00 10. 8 225 

Max Energy Flux 

(q ,) (kW/cm 2 ) 

raa 

max 

1.9 3.6 11.0 

Av. Energy Flux 

q , (kW/cm^) 
rad 

1.19 2.31 7.16 

% Power Loss 

■H 
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in the dissociated and ionized state. We have therefore made calculations 
including these properties, which play an important role in determining 
the heating, since the temperature variation, and therefore the property 
variation, across the boundary layer is so large. 


A simple set of boundary layer equations is obtained in the standard 
way by transforming from the wall - oriented coordinates s, y along and 
normal to the wall to the similarity coordinates , T| . This transforma- 
tion is given in terms of velocity u g , viscosity}! , and density P e at the edge 
of the boundary layer by 


/. (a) = / P V u ds, H (s,y ) = % - 

*; JIT 



(8.4) 


The dependent variables velocity u and total enthalpy H are written in 
normalized form as 


u/u e = f^ , H/H e = g (8.5) 

where the subscript 'H denotes differentiation with respect to Tl , The 
boundary layer momentum and total enthalpy equations are transformed 
using (8. 4) and. (8* 5). The "local similarity 11 assumption is that 1 and g 
depend only on T ; the equations are then 

(8* 6a) 


(8* ob) 


( 8 . 7 ) 

and k denotes the thermal conductivity, |! the viscosity. 

The hypothesis of "local similarity" has enabled us to convert the 
partial differential equations of momentum conservation parallel to the wall, 


(C f ) 4 f I 4 

iw Tvn 


* (°e z\ 

P [— -l = 0 
\ a r\ J 


*1 e 


2 - 

h~“ TtT [ C (* " ^ ^ f *n-n j =0 


where we have defined 


Pr = c u/k , C = 0}!/o u , R = 2d In u /din I 
p e e e 
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and energy conservation, into ordinary differential equations in the variable 
T) . The variable £ , which depends only on distance along the wall, s, by 
(8.4), appears as a parameter in u and in (3*, which involves the derivative 
&( tt e with respect to £ . The gas property functions Pr, C and 0 /p depend 
only on the state of the gas, so are functions of p and T (or p and enthalpy 
h = H - u e ^ 1^/2) and so can be written as functions of g at each value of 
£(s). Thus, when the properties external to the boundary layer are speci- 
fied as functions of s, Eq. (8. 6) can be integrated to find f(Tji), g(TV), 

Two further simplifi cations are useful. Since u e /H e is small, and 
Pr is near unity for hydrogen over the range of interest, we will put Pr = 1 
in (8* 6b). Further, the third term in (8. 6a), the pressure gradient term, 
is small so we will approximate @ e /p by h/h , When this is expressed in 
terms of g, a new pressure gradient parameter $ appears, 


H H d In u 

R = S. ft* ~ ? — g 

g h 8 2 b din l 

e e 


Now the boundary layer equations are 


“t,n + S< ®* V> 5 8 


( 8 . 8 ) 


(8.9a) 


<Cg > + fg - 0 (8.9b) 

T1 ' 

The boundary conditions are no slip at the wall, a given wall 
entttalpy h based on the wall temperature T^, and approach to the exter- 
nal conditions at large Tl : 


^ (0) = ^ S (0) = h w / M e = 6 W ’ f ri ( “ } = g (m) = 1 < 8 * 10a > 

For the case where there is no flow through the wall, we have the fifth 
condition that the normal velocity is zero there, which is given by 

f (0) = 0 (8.10b) 
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When the external conditions p e , ta e , \i & are given from the core 
streamtube solutions, and the p and p, are defined in terms of h, then 
(8. 9) and (8. 10) can be integrated. The p (h) function was found from the 
equilibrium properties for hydrogen given in Ref* 33, as was the function 
h(T), The viscosity was found from Ref, 3?, as a function of T and p, 
completing the necessary information. The integration was performed by 
the method of quasi -linearization. 


The results of the integration for T e = 14, 000 K p = 3 atm, are 
shown in Fig , 8. 3 as plots of the shear and heat transfer parameters 


f^(Q) and (0) against the pressure gradient parameter p, for various 
values of g w . These are related, to the shear stress t w and heat flux q w 
by 


T - 
W 


2 

P |! u ' 
w we 


/rr 



0 |i u H 
w w e e 

J Z% Pr 




( 8 . 11 ) 


In addition to the solid curves for equilibrium hydrogen properties, 
the dashed curves show the results of assuming that the p p product does 
not vary across the layer, i e, , C = 1. For the small values of wall to 
external enthalpy ratio g w in which we are interested (g w « 0. 02 for 
T w = 1, 300 K, the melting point of copper) it is very important to consider 
the variation of C, This is not surprising, since for 3 atm, T e = 14, 0'0OK f 
T w = 1,30'OK, C varies from unity (by definition) at the outer edge of the 
boundary layer, to 7,4 at the wall. Another notable feature is the very 
slight dependence of (0) on pressure gradient parameter (3, especially 
for g. w — 0*02. 

Based on calculations like those shown in Fig* 8.3, and the use of 
Eq. ( 8. 11), the heat transfer distribution along a 1, 30'0K wall in contact 
with the plasma core was calculated for the 3 atm, 10 kW ease. Because 
conditions in the core for the 10'0 kW and 5 000 kW eases are only slightly 
different, the resulting heating rate is also a good estimate for these two 
cases at 3 atm. The distribution of q w is shown in Fig. 8.4, and gives 
fluxes from 10 kW/cm^ to a peak of 34 kW/crn^ near the throat. If these 
values are integrated over the surface, the total convective heat loss rate 
is found from 


Q = / 2tr R q dx (8,12) 

c */ w ; 

G 

The values that result for the 10 kW, 100 kW and 5000 kW cases (each of 
which have a different radial size), for integration from x - 1. 5 to x^, are 
17,53 and 374 kW respectively. 
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It is immediately evident that this uncoupled calculation produce 
unrealistic results for both the IQ kW and 100 kW eases, since the calculated 
heat loss rate is a large fraction of the total power. Such a loss rate would 
decrease the core temperature, which in turn would decrease the loss. 

For the low power eases, then, the convective loss must be coupled to the 
core calculation for a realistic model of the fLuid mechanics. For the 
5 00 kW case, however, the loss calculated is only 7. 5% of the laser power, 
which is probably not a large enough loss to require coupling to the core. 

For this case, the heating distribution shown in Fig. 8.4 is probably realistic. 
This certainly indicates t'he need for some protection of the rocket walls, 
since rates of 30 kW/cm.2 cannot be handled by simple means. However, 
the larger rates do occur near the throat, where the stream-tube is of small 
radius, so the area to be protected is not large. The throat radius for the 
3 atm, 5000 kW case is about 1 cm. 

Realistic heat fluxes for the lower power engines can only be 
determined by a coupled calculation of the core and boundary layer, which 
was not undertaken under the present contract. 

The boundary layers for the 10 atm cases were not calculated, but 
for Laminar heating, the heat fluxes should scale like the square root of 
pressure, increasing them by a factor of about 2. But the areas of the 
streamtube surface are nearly 10 times smaller because of the short 
length, so the values of Q G should be about one fifth of those for 3 atm. 

This means the coupled calculation is still necessary for the 10 kW case 
at 10 atm, might not be necessary for the 100 kW case, and of course ljj 
not necessary for the 5000 kW case. The heat protection of the walls 
will be even more important though, since q^ is higher by a factor of 2, 

One method of heat protection is the injection of cold gas normal to 
the wall, which alters the temperature profile and blocks some of the 
heat. Some sample calculations have been made to illustrate the effect 
of this method of reducing the heat flux. The only change in the calcula- 
tion is the Introduction of a non-zero value of f (0) to replace the boundary 
condition (8. 10b), because now the mass flux through the wall is not zero, 
but is related to f (0) by 


(pv) w = -p w fci w f (0) / 72 5 


(S.13) 


This introduces another parameter into the equations, f (0). The results 
of two calculations, at stations x = 1 and x = 2. 7, are given in Fig. 8, 5 as 
curves of g (O') against f (0), The left end of these curves is the zero 
injection or 1 ^ olid wall case f (0) = 0. It can be seen that the heat transfer 
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parameter (0) drops sharply as £ (0) decreases, which mea^s the 
injection velocity increases by (8* 13). If we look at injection with 
f (0) = - l t g^ (0) at x = 1 is reduced by a factor of 4, leading to a heat 
flux of about 2*5 kW/cm^ for an injection rate of 0.014 gm/cm^-s, 
from (8* 13). At x = 2, 7, g^ {0) is reduced by a factor of 2, yielding 
17 kW/cm^ with an injection rate of 0.017 g/cm^-s. 

To see the magnitude of this rate of injection, one may compare 
it to the mass flow In the core, which is 0, 3316 g/cm^-s over the initial 
area. Thus for the 10 kW ease, the core mass flow is 9. OE-3 g/s, while 
an injection of 0.015 g/cm^-s over the surface area of the streamtube 
up to the throat is 22E-3 g/s, more than twice the core mass fLow. On 
the other hand, for the 5000 kW ease, the cor e mass flow is 500 times 
larger, 4.5 g/s, while the area is only A /500 = 22 times larger, so the 
injected mass would be about 0.49 g/s, or only 11% of the core mass 
flow. The surface to volume effect favors injection in the large engines, 
and it appears to be one possibility for reducing the wall heat load by a 
factor of 2 or more. Notice also, that large reductions in heat flux are 
only needed in the limited regions of high heating, not over the whole 
streamtube surface, so injection is not needed over the whole nozzle 
surface, which reduces the mass which must be injected. 

The purpose of the present calculation was to show how much 
heating reduction can be achieved by mass injection. Considerable addi- 
tional calculation would be needed for an engineering analysis of the 
proper amount and distribution of injection for efficient heat protection. 

8. 4 ABSORPTION BY PARTICLES IN BUFFER GAS 

It appears that for specific device configurations the radiation 
escaping from the plasma will be of sufficiently high intensity that some 
sort of shielding mechanism will be required in order to maintain wall 
integrity. One possibility would be to introduce small carbon particles 
into a buffer gas flowing around the plasma. In this event the carbon 
particles would intercept the plasma radiation, heat up and conduct heat 
energy to the surrounding buffer gas. These particles would also radiate 
with an intensity corresponding to their own characteristic temperature 
and could further reach vaporization temperatures depending on the 
particle flow time and their radiative flux loading. 

In the following analysis first order estimates are provided for 
the particle mass loading required, to protect the walls against plasma 
radiation. The particles are assumed to have constant absorption proper- 
ties, velocities and size. An understanding of the phenomenology of the 
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plasma/particle/buffer gas interaction, including the all-important 
question of particle survival, will require a significantly more 
sophisticated analysis of the mixing region between the plasma and 
buffer gas flows. When such an analysis is performed the required 
particle mass loadings can be more accurately defined. 


Small carbon particles, of the order of 1 p.m, are very efficient 
absorbers of radiation in the visible and near infrared wavelength 
regions, approximately acting as black-body absorbers. A large fraction 
of the hydrogen plasma radiation will unfortunately originate in the far UV 
where little is known about the absorptive properties of carbon. In the 
following analysis it will be assumed that the carbon particles absorb as 
black bodies at all wavelengths. 

The configuration considered is a cylindrical plasma of radius r ^ 
surrounded by a particle seeded buffer flow of thickness /\r. The 
particles are taken to be spherical with radius Rp and to have a constant 
number density np in the buffer flow. Each particle has an absorption 
cross section of ttR ^ and thus the optical depth, t, across the buffer 
layer is defined by ** 


T = n frR ^ Ar (8. 14) 

p p 

An optical depth of 3 Is sufficient to absorb 95% of the plasma radiation 
and this condition is met when the particle number density is given by 


n 




s* ** 


(8.15) 


It is perhaps more interesting to examine this quantity in terms 

of the total mass /cc of particles in the buffer flow, M , 

F 


M =n p (4/3) ttR ' 
P P P P 


(8. 16) 


where p is the specific density of carbon which is ~2 gm/ce, 
Eq. (8. 1*1) and (8. 16) results in 


R 


M p = 4 p p 


hr 


R 


A* 


(8. 17) 
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Equation (8. 17) ean readily be extended to deduce the mass flow 
of particles required to absorb the plasma radiation. 


where 


M = M u, A, 
p p b b 


A fe = rr (2 r 1 Ar + (A r) ) 


{8. 18) 


(8. 19) 


and is the buffer flow velocity. Since it is desirable to maintain a con- 
figuration such that A r / r i « 1# the approximation 


A b ~ 2 nr 1 Ar 


( 8 . 20 ) 


is made. Combining Eq. (8.17), (8. 18) and (8.20) results in 



8TTp p U b R p r l <s48 


u, R r , 
b p 1 


( 8 . 21 ) 


It can be seen from Eq. (8. 21) that as long as fir « rj, the 
particle mass flow required to absorb the plasma radiation is Independent 
of the buffer flow thickness and linearly proportional to both the particle 
and plasma radii. In the interest of optimizing thrust one wishes to 
minimize the buffer flow; hence, it is clear that the smallest particle 
size consistent with the assumption of biackbody absorption should be used. 

Another interesting feature is that since the required buffer How 
is Linearly proportional to plasma radius, the ratio of buffer flow to plasma 
flow is size dependent; i. e. 


M 8 p R u. 

P _ H P P b 


M 


pi 


P , r. u , 

V i pi 


( 8 . 22 ) 


where u ^ is the plasma flow velocity. For example, for a 14, Q00°K hydrogen 
plasma M a pressure of 3 atm. 



10 


_ R u, 
7 P *> 


r. u . 
1 pi 


(8.23) 
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Note that relationship (8.23) implies that for small devices, L. e. 

< 1 cm, the buffer mass loading would of necessity be orders of 
magnitude larger than that of the plasma for reasonable particle sizes 
unless ub « Thus these preliminary numbers suggest that particulate 

absorption may not represent a viable technique for wall protection in the 
flowing plasma configuration. Note that even if the buffer flow had a 
velocity two orders of magnitude lower than the plasma, a condition 
difficult to maintain along the flow length, there may still be an unacceptable 
constraint on particle size. 

Although the analysis presented above is admittedly crude it would 
suggest that particle shielding wouLd provide too Large a weight penalty for 
proposed operating conditions. A more detailed analysis would be required 
to confirm this. Furthermore the possibility of operating under different 
conditions, for instance at higher plasma pressures, should also be 
cons ider ed. 
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9. PERFORMANCE OF FLOWING CORE ROCKET 


The performance of the flowing -core Laser -heated rocket is affected 
by many more variables than the ones considered in this report. Most of the 
effort reported here has dealt with the fluid mechanics, laser radiation 
absorption and plasma emission of the flowing core, and the convective heat 
transfer from the core stream tube. Based on this work, only the performance 
obtained from the core stream tube can be obtained, that Is, the thrust and 
specific impulse of the core gas when it is expanded to some exit condition. 

The resulting values will be high estimates because no Losses past the throat 
are considered, nor has the convective heating loss been properly coupled to 
the core streamtube for the lower power Levels as yet* On the ocher hand, 
added thrust which might be supplied by a buffer gas Layer used for heat 
protection is not included either. With the understanding of the partial nature 
of the performance which can be calculated here, we can describe the method 
for the calculations. 


We begin at the throat of the streamtube with the enthalpy h^ and 
velocity u^ found from the program described in the section on the Flowing 
Core Streamtube. The flow is then expanded isentropically to a specified 
exhaust pressure p ex by linear interpolation in hydrogen property tables in 
Ref. 33. This provides the exhaust conditions of enthalpy and temperature 
h ex , T ex , Conservation of total enthalpy during the expansion (no losses) 
then provides the exhaust velocity as 


u 

ex 



- h ) 
ex 



The thrust is then 


Th = <p 


ex 


p ) A + 
am ex 


m u 


ex 


(9.1) 


where P am is the ambient pressure. We will give the so-called ideal thrust, 
for which p ex = p arn * so that 


Th = mu , I 

ex sp 


ex 


(9.2) 


give the thrust and specific impulse, where g is the acceleration of gravity. 
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The conditions at the throat, as calculated by the core streamtube 
program, the mass flow rate from Table 7,2 and the entropy s from 
Ref, 33 are given in Table 9* 1* The exhaust conditions, I S p and Th for 
Pex “ 1 atm, as found from Ref, 33, and Eqs, (9. 1) and (9.2) are given 
in Table 9.2, The maximum thrust is obtained by expanding to p ex - 0, 
since all energy is then converted into velocity, h ex being zero. This 
limiting ease is given in Table 9.3, as obtained from Eqs. (9, 1) and (9.2) 
by setting h ex = 0, 

Table 9*2 shows that the expansion to 1 atm of the 3 atm cases is 
really a very small expansion, since the throat is at about 1. 7 atm. There- 
fore, the gas in the exhaust is still very hot (near 13, 000 K) and its 
enthalpy is only a few percent less than at the throat. It is this few per- 
cent which is converted to additional velocity, so the increase from the 
throat velocity to the exhaust velocity is also only a few percent. However, 
this velocity is sufficient to produce I S p in the range of 1400 to 15 00 s. 
Because of the small mass flow rate, the thrust is in the range 0*1 to 6b N 
(0, 03 to 15 libs , ) , 

The 10 atm eases have a greater expansion ratio, since the throat 
is at 5, 6 to 5, § atm, so they have a greater conversion of enthalpy to 
velocity, exhausting somewhat cooler gas, at speeds nearly double the 
throat speed. They then produce higher values of I , Ln the range 2400 
to 2500 s, and somewhat larger thrusts, 0,2 to 100 N (0,05 to 25 lbs. ), 

For the maximum possible expansion, in Table 9. 3, we see much, 
less difference between the 3 atm and 10 atm cases, since the expansion 
ratio is infinite for both* They all yield I S p in the range 3000 to 3300 s. 

In fact, the 100 and 5000 kW cases have nearly the same at I g p at both 
pressures. Only the highest percent loss case, 10 kW, has noticeable 
difference in I S p with pressure. The thrusts for the vacuum exhaust 
cases are also higher than the atmospheric exhaust cases since thrust 
is proportional to u ex . 

Even though a number of loss mechanisms were omitted from 
this calculation, and they will reduce the performance calculated here, it 
seems quite clear that values of I S p above 1000 s can be obtained, even, 
over a 3 atm to 1 atm expansion. 

The thrusts produced may seem rather small, because of the very 
small mass flow rates. However, the Th/P ratio for the vacuum expansion 
core is in the range 20 to 30 mN/kW, which seems to be typical for electric 
propulsion schemes, to which the laser -heated rocket is the most comparable. 
Of course, the laser -heated rocket does not have to carry its own power 
supply, so its weight for a given mission should be smaller than an electric 
propulsion engine. 


- 106 - 


- 401 - 


TABLE 9. 1 

THROAT CONDITIONS 


P(kW) 

p (atm) 
o 

h th( J / kg ) 

%‘ N ^ 2 * 

T th« K) 

s th (J/kgK) 

U th (m/S) 

• , 2 
m(kg/m 

10 

3 

7. 99E8 

1. 73 E 5 

13490 

2. 1721E5 

1.25E4 

9.02E-4 

5000 

3 

9. 63 E 8 

1. 71E5 

14280 

2. 2940E5 

1. 33E4 

4. 51E-1 

100 

3 

8. 97E8 

1.73E5 

13990 

2. 2456E5 

I.29E4 

9. 02E-3 

10 

10 

8. 99E8 

5. 67E5 

14970 

2. O900E5 

1. 34E4 

9. 02E-4 

100 

10 

9. 46E8 

5. 77E5 

15240 

2. 1485E5 

1. 34E4 

9. 02E-3 

5000 

10 

9. 73 E 8 

5. 82 E 5 

15380 

2. 1646E5 

1. 34E4 

4. 51E-1 
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TABLE 9. 2 


PERFORMAi CE FOR p atm 


P(’.W) 

(atm) 

h ex ( J/kg) 

T (K 
ex 

10 

3 

7. 79E8 

1 3000 

50:00 

3 

9.40E8 

13760 

100 

3 

8. 74 E 8 

13470 

10 

10 

6. 73 E 8 

12350 

10:0 

10 

7.47E8 

1 2800 

5009 

10 

7. 67E8 

12930 


ex 


m (m/s) 
ex 

I Js) 
sp 

Th( N) 

Th{lb) 

1. 38E4 

1410 

0. 124 

0. 0279 

1.47E4 

1500 

66. 4 

14. 9 

1.44E4 

1470 

1. 30 

0. 292 

2. 48E4 

2530 

0. 224 

0. 0504 

2. 37E4 

2420 

2. 14 

. 482 

2. 40E4 

2450 

108. 

24. 3 






~ i 
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TABLE 9. 3 


t 

f 



P(kW) 

P o (atm) 

u (m/s) 
ex 

I (s) 

Sp 

Th(N) 

Th(lb) 

10 

3 

2. 96E4 

3020 

0.266 

0. 0599 

50:00 

3 

3. 24E4 

3300 

146. 

32. 5 

i! 00 

3 

3. 13E4 

319© 

2. 82 

0. 634 

10 

10 

a. 14E4 

3200 

0. 283 

0.0637 

100 

10 

3. 22 E4 

328'0 

2. 90 

0. 652 

5000 

10 

3 . 26E4 

3330 

147. 

33. 1 


I 
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10. STATIONARY CORE ENGINE 
10.1 INTRODUCTION 

An alternative to the flowing plasma core type of rocket is one 
borrowed from the concept of a gas -core nuclear rocket. In that type 
of device a gaseous core of fissionable material is maintained in a 
stationary region in the center of the rocket while fissioning at high 
temperatures. It is surrounded by a flowing propellant gas which confines 
the reacting core, and is heated by thermal radiation and convective 
processes. The nuclear fuel hums above 100, 000 K, and it must be kept 
at least at atmospheric density to maintain criticality, so it operates at a 
core pressure of several hundred atmospheres. 

An analogous laser -heated rocket would have a stationary plasma 
core, which absorbs the laser radiation and is thereby maintained at high 
temperature. In turn, this core re -radiates the absorbed laser energy 
to the surrounding propellant gas, which serves to confine the core and 
produce thrust, just as in the gas -core nuclear rocket. This concept may 
be compared with both the gas -core nuclear rocket and the flowing core Laser 
heated rocket. 

Compared to the former, there is no criticality requirement 
imposed on the core. However, Laser energy absorption, which is the 
mechanism for putting energy into the rocket, is dependent on the core 
density. The more dense the core, the better it absorbs, and the shorter 
it can be in the direction, of the laser beam. But there is more freedom in 
the range of densities, and therefore pressures, in the core for the laser - 
heated case than for the gas -core case. Another difference is in the 
amount of core plasma which can be allowed to flow out the rocket nozzle. 
For nuclear rockets, one would like to have as little core material as 
possible entrained by the propellant gas and Carried out the nozzle, since 
it represents a loss of energy of fission as well as of radio-active material. 
In the laser rocket, "detrainment" of core material is not so serious, 
although it does represent a loss of propulsive efficiency, since it will 
leave the rocket at rather high temperatures thus wasting energy. However, 
if we let too large a fraction of the core flow out with the propellant gas, we 
are back to the flowing core concept. Preventing this detainment is not so 
easy in the laser rocket, since the core may not be of high density relative 
to the propellant gas, as in the nuclear core. 
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Compared to the flowing -core laser -heated rocket, there are 
several new physical mechanisms which enter* One is the detrainment 
problem referred to above, which acts against the confinement of the core* 
The densities of the core and propellant may be nearly equal, which makes 
confinement difficult. Further, the heating of the propellant and its sub- 
sequent density decrease occur while it is flowing past the core, making 
detrainment easier as it flows. A second important mechanism is the 
efficient emission of radiation by the core, and its efficient absorption by 
the cooler propellant gas. In the flowong -core scheme, emission by the 
core is to be minimized as a loss of energy from the primary propellant, 
while in the stationary -core this emission is the main process by which 
energy is input to the propellant. Likewise, absorption by the propellant 
gas Ls now a primary source of propulsion rather than a necessary heat 
protection device, as in the flowing -core rocket. In fact, one of the main 
complications of the stationary-core rocket is in the control of the flowing 
propellant gas. This is achieved primarily by control of the flow area 
transverse to its direction of motion, i. e, , the cross -section area of 
the propellant gas streamtube. The same Lb true, of course, in the flowing- 
core case; its flow properties are largely controlled by the cros s -sectional 
area. But there Ls one crucial difference. In the flowing -core ease, the 
propellant gas absorbs energy from a laser beam propagating in the axial 
direction, so its absorption properties depend on its axial size. Thus the 
control of flow and absorption depend largely on two independent geometrical 
parameters, the cross -section and the axial length, However, in the 
stationary -core rocket, the propellant gas absorbs its radiant energy 
mainly from emission in the transverse direction, the same direction which 
controls the flow. This means that flow control will be intimately connected 
with absorption, and the propellant gas flow area must be designed with 
regard for both the axial pressure, velocity and density distributions and 
for the gas depth for efficient absorption. 

It appears, then, that the analysis of the stationary -core rocket is 
considerably more complex than that of the flowing -core rocket, with 
increased emphasis on plasma emission and subsequent absorption, and 
a more complicated interaction between flow and radiation. 

10,2 OSDEft OF MAGNITUDE ANALYSIS 

The quantitative evaluation of a stationary-core rocket can be 
started by using some simple relations among the geometric and physical 
parameters which serve to indicate the physics which goes into such an 
analysis, as well as obtain some order -of -magnitude estimates of the 
parameters . 
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Taking a cylindrical core of radius R c and length L, the laser 
power P and intensity I are related by 


P = ttR 2 I (10. 1) 

G 

The absorption of laser energy occurs with absorption coefficient k^,, a 
function of temperature T. The volume emission by the plasma, P c , 
also depends or. T. In steady state the plasma must emit the full power 
P, so that 


P = ttR LP, 


( 10 . 2 ) 


Equations (10.1) and (10.2) give 


L = I/P 

c 


(10. 3) 


The absorption and emission in the core are related by the steady state 
requirement that the absorption k^I equals the emission F e : 

k I = P (10. 4) 

J— i G 

Equations (10, 4) and {10, 3) show that with these estimates, the core 
length L is equal to the absorption coefficient k . 

w * 1_# 

Of the emitted radiation per unit volume P e some portion P G <P e 
wiLL be absorbed by the propellant gas , depending on the depth of this gas, 
and its absorption coefficient kp. If nn^ is the mass flow rate of the 
propeLlant gas , the energy per unit mass added is rr R e ^ P*,/ rinp . If its 
initial enthalpy is hi, then its heated total enthalpy is 

H = h. 4 uf/2 + ttR 2 Lf^/ra (10.5) 

h i 1 c e P 

This gas is then expanded at constant entropy to the chosen exhan-st pressure, 

p , yielding an exhaust enthalpy h , and an exhaust velocity 
fix ex 


U ex J 2 (H - h ) 
h ex 


(10.6) 


and associated specific impulse 


*ex/g = 1 


sp‘ 
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With these relations, estimates of the parameters of a rocket 
can he made. Suppose we give the power P, core temperature T c , core 
pressure p B , exhaust pressure Pex and specific impulse desired, which 
defines uex« From T G and the radiation properties of hydrogen we can 
find P e and kj^. Then (10.4) gives I, (10,3) give L and (10. 1) give Rc» 
so the core size is determined. 

The mass flow rate of propellant m p may now be found from u ex . 
Equation (10.6) gives Hh - hex* H we assume uj^/s is a small contribution 
to this enthalpy difference, then we know the desired difference between 
the enthalpy of the propellant gas at the end of the heating zone, which we 
will take to be at constant pressure p e , and its enthalpy at the exhaust 
pressure P ex . These two states are connected by constant entropy. The 
conditions uniquely define both states, so the values of T^, hfo, T ex , h ex 
can be found. Again assuming the velocity contribution to the total enthalpy 
is small, we take hh - h^ to be energy added to the gas. Now we estimate 
the volumetric absorption P^ as a fraction of the emission P c . Then 
(IQ. 5) yields the propellant mass flow as 

m p = rrR e Z L F^/d^.h.) (10.7) 

where the initial enthalpy is assumed known. The ideal thrust can also be 
found now as 


Th = m u 

P ex 


( 10 . 8 ) 


To cheek the neglect of the velocity compared to the enthalpy during 
heating, and to estimate detainment effects, we need to define the area in 
which the propellant gas flows. This is in turn related to the absorption 
length of plasma radiation in the gas, and so depends on the detailed proper- 
ties of the gas. If the radius of the outer edge of the propellant gas is R 
at any station, its velocity is 


u p = m p /p p rr(Rp -R e 2 ) < 10 ‘ 

and p p is known from the given gas state at the inlet, and the end of heating 
state, which was found above. 

The detrainment occurs in the shear layer between the stationary core 
and the flowing propellant, where the two gases mix. If the thickness of the 
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mixing layer on the propellant gas side is 6, then the core gas mass flow 
there is 


■ 

m 

e 



a' dy 


( 10 . 10 ) 


I I 

where p , u are the density and velocity of the core gas in the mixing 
layer. By normalizing, this can be written 


. r 1 

m. = 2 TT R $ P b_ I 
c e c P J 

0 

The factor in front gives the order of magnitude of the mass flow rate out 
of the core at any station, and the total core de trainmen t. 

The thickness 6 Is determined by a balance between the diffusion 
of core gas into propellant gas and the convection of the propellant pa st th e 
core. This can be estimated by combining the static diffusion length t 
with the time it takes for the propellant gas to flow past the core, 

t - ^*dx/up. Here D is the diffusion coefficient of core gas into 
propellant gas. Thus at the end of the core 


P' 

P 


u' 

U. 




( 10 . 11 ) 


5 = 



(10. 12) 


To evaluate the integral, we will use a linear increase of up from u-p. to 
u Ph> so 


6 


/ DL 
V u^, -u_, 

\ Ph Pi 


In 


a Ph\ 1/2 
U Pi / 


(10. 13) 


Now the factor in fr ont of the integral in Eq.. (10. 11) can be calculated as 
an estimate of detrainment. 


10.3 NUMERICAL EXAMPLE 


As an example of the use of these formulas for estimating the 
performance of a small stationary core engine we may take as a core: 


P = 10 kW, p. = p = 3 atm, T = 14, 000 and 18, D00K 
1 r c c 

p ex = 1 atm, I sp = 1000 s 

Then the core conditions are: 


t g (K) 

14, 000 

18, 000 

p„ (kW/em 3 ) 

56.4 

200 

(erti' 1 ) 

1 

3 

I (kW/em 2 ) 

56 

67 

R e (cm) 

.24 

.22 

L (cm) 

1 

. 33 

P c (g/em 3 ) 

2.33 E-6 

1.25 E-6 

p v (kW/em 3 )' 
0 

29 

90 

frR| L F* (kW) 

5.2 

4. 5 


The emitted radiation which is assumed to be absorbed by the 
propellant gas, P* , is that below 0, 1 |i-m. The remaining energy is lost. 
This is a rough estimate based on our present knowledge of absorption 
by hydrogen, and ignores IR absorption. 

Turning to the propellant gas, I fi p = 100*0 s corresponds to 
u ex - 980 E5 cm/s and a kinetic energy per unit mass of 4. 8 E 7 J /kg. 
Hydrogen properties at p ex = 1 atm then show that at T ex = 5000 K, 
hex r 0* 93 E8 J/kg at an entropy of s ex = 170 J/kgK, which leads to 
Hex - 1 • 41 e* J/kg, The corresponding point at p e = 3 atm is at 7000 K , 
t% = 1. 43 E8 J /kg, = 170 J/kgK, pj^ = 5. 3 E -6 G g/em^, which gives 
the state of the propellant gas after it is heated. Knowing the power to 
be absorbed to heat from = 300 K, = 2. 10 E8 J/kg, P^ = 3 atm to this 
state, we find the mass flow rate of propellant gas, mp, as 0,015 and 
0,013 g/s for the two core temperatures. 

To find the propellant gas speeds and core detra Lament we must 
now make some assumption about absorption depths to absorb the UV radia- 
tion A reasonable estimate Is 0. 1 cm, leading to Rp = R G + 0. 1. 
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For this constant area flow of propellant gas, the value of up p p comes 
from Eq. (10. 9), and since the density is known from the constant 
pressure and known temperatures Ti “ 30GK, - 7OO0K, the velocities 

can be found, as can the constant factor in Eq. (10, 11), using (10* 13) 
for 6. The diffusion coefficient D is evaluated at = 7000 K by scaling 
the value for H, diffusion in Reference 32 down to 3 atm. 


T c (K) 

14,000 

18,000 

wip(g/s) 

l.SE-Z 

1.3E-2 

Pf> ip(g/em 2 -s) 

8 , 2E-2 

7.7E-2 

u pi (em/s) 

332 

312 

u Ph (cm/s > 

15. SES 

I4.5E3 

D(em 2 / s) 

230 

230 

§ (cm) 

0.23 

0.14 

21TR 6p u (g/s) 
c c c JF n 

1.2E-2 

3.5E-3 

Core Mass (g) 

2.9E-7 

6.3E-3 

(Core ma8s)/m (a) 

2.4E-5 

1 .8E-5 


The first observation from these results is that the kinetic energy 
in Up., i E4 j/kg is small compared to h^ -1*43 E8 j/kg, and so was 
indeed ignorable. The next observation is the tremendous acceleration 
of the propellant gas while being heated, although it does not reach the 
sonic speed, which is a^ = 9. 5E5 cm/s, but stays at very low Mach 
number* Finally, we see that the detrained core mass flow, by our 
estimate, could be comparable to the propellant mass flow. Further, it 

C 

empties the core in a time of order 10~ 3 seconds. While this estimate of 
detrainmenfc is quite crude, and may be on the high side since the diffusion 
coefficient was evaluated at Tft, still it shows that detrainment r an be a 
big potential problem for the stationary-core concept* 

This example calculation is meant to illustrate the sort of analysis 
necessary to define the properties of a stationary -core laser -heated 
rocket. Most of the uncertainty in the results comes from Lack of 
detailed knowledge of the absorption proees ses in the propellant, and of 
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the core detrainment process. Both of these areas need further study 
if the stationary -eo re concept is to be pursued. With some more detailed 
models of these two processes available, parametric studies could be 
done at various values of T e , p G , P. There may well be preferred 
regions in which to operate such an engine, and parametric studies 
might uncover them* 


11. CONCLUSIONS 


Based on the studies of laser-heated rockets described in this 
report, as well as the knowledge gained in the previous study reported 
in Ref. 1, we can draw some conclusions about the feasibility and opera- 
tion of these devices. We can also identify important areas where our 
knowledge is incomplete, and further research is needed to put the con- 
clusions on a sound scientific basis. 

1. The flowing core rocket starts with a laser supported 
combustion (LSC) wave. For a given laser power, the cross-sectional 
area of the rocket engine is inversely proportional to the laser intensity 
at which this wave operates. The temperature level of the pLasma core 
and the mass flow rate through the engine are also determined by the 
LSC wave, so a knowledge of its properties are vital to the definition of 
the rocket. However, there is at present no theoretical or experimental 
information on such waves in hydrogen, and even in air they are not 
well-understood. The properties used here were estimates based on 
extrapolating a simple air theory to hydrogen. However, LSC wave 
properties in hydrogen represent the major uncertainty in the present 
study, and should be studied both experimentally and theoretically to 
increase confidence in the feasibility of laser-supported rockets using 
hydrogen. 

2. Performance of the flowing -core streamtube as a rocket is 
good. Values of I S p for expansion to l atm exhaust are 1500 to 2500s 
while expansion to a vacuum leads to 3000 to 3300 s. 

3. Only a limited number of calculations were made, using 
three laser powers at each of two pressures. The LSC wave properties 
and the core velocity distribution were not varied. To map out the region 
of usefulness of laser-heated rockets, more extensive parametric calcula 
tions are needed, 

4. Breakdown of hydrogen to produce an electron density 
which is comparable to that of thermal ionization can be produced by 
either an electric are discharge or by a focused laser beam. Whether 
this will actually initiate an LSC wave has not been determined, and 
probably needs to be verified by experimental investigation. 
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5. The radiation emission of hydrogen has been put into a 
form useful for fluid dynamic calculations. The black body portion of 
the emission depends crucially on the radial temperature profiles, and 
has not been included in the present one -dimensional calculations. Its 
inclusion would require some model for the radial-variation of temperature. 
It couLd be a s LgnLficant effect, however, and should be investigated 
further. 

6. Based on the transparent radiation emission, radiated energy 
loss is not a substantial fraction of the laser energy input, varying from 
20% for smaller engines to 4% for larger ones. This conclusion Is 
sensitive to the operating temperature of the plasma core, and could 
change if this temperature is even 1000 K higher than calculated here. 

7. The energy flux radiated is sufficiently high to require some 
kind of active heat protection for the walls, ranging up to 75 kW/em^. 

8. The convective energy losses calculated by assuming a 
hydrogen boundary layer at the edge of the plasma core, without coupling 
the Loss to the core flow, are too high to be realistic for aLl but the 5000 IcW 
laser powe r. A coupled calculation, not undertaken here, is needed for 
smaller power levels* 

9. Far the 5 000 kW power, the convective heat losses are of the 
order of 30 kW / em^ near the throat, at 3 atm, and perhaps twice that for 
10 atm, so they add significantly to the need for heat protection, though 
they only occur in a small region. 


10. Normal injection of hydrogen into the core streamtube seems 
to be a useful way of reducing the convective heat losses, based on a few 
sample calculations. Further calculations are needed to establish this as 

a practical method of heat protection. Tangential injection was not studied, 
but should also be investigated for heat protection. 

11. First estimates indicate that particle absorption does not 
seem a useful way to protect the walls from radiation, because of the high 
mass loading of particles required. This subject needs further study. 

12. The length of the plasma core, up to the throat, is controlled 
by the inverse Bremsstrahlung absorption coefficient. Larger absorption, 
caused by larger pressures, leads to shorter lengths. At 3 atm, this 
length is 4 cm, while at 10 atm it is 0.5 cm* These short lengths it the 
high pressure lead to short, stubby streamtube shapes, which may be 
undesirable from a fluid -dynamical point of view, especially for larger 
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powers. This imposes a Limitation on how high a pressure level can 
be used for operation. 

13. The stationary- core rocket engine concept was studied 
qualitatively and semi -quantitatively. The physical processes involved 
in this concept are more complex than those of the flowing-core concept. 
Two major uncertainties were uncovered. One is a lack of knowledge 
of absorption lengths for high temperature hydrogen radiation in low 
temperature hydrogen. The other was a lack of knowledge of how the 
flowing propellant gas entrai ns core gas. Further study of these two 
areas would enable estimates of stationary-core rccket performance to 
be made with more confidence. 
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